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I.  SLIN^MARY 


A.  PURPOSE  AND  SCOPE 

In  this  report  we  discuss  the  problem  of  designing  optimum  signals 
and  receivers  for  an  active  sonar  system  operating  in  an  environment  in  which 
the  principal  cause  of  interference  is  reverberation.  We  have  considered  only 
the  problem  of  target  detection  (as  opposed  to  the  problem  of  range  or  Doppler 
estimation) . 

In  many  sonar  systems  the  primary  limitation  on  performance  is 
energy  scattered  back  to  the  receiver  from  various  objects  in  the  ocean.  If  the 
scattering  structure  were  known  in  detail,  the  reverberation  could  be  treated  as 
a  deterministic  return  and  the  signal  and  receiver  designed  to  eliminate  it.  Be¬ 
cause  of  variations  in  the  ocean,  however,  this  type  of  model  quickly  becomes 
unrealistic  or  unmanageable.  For  this  reason,  we  chose  an  alternate  approach 
in  which  the  reverberation  return  is  treated  as  a  random  process.  Characteriz¬ 
ing  the  reverberation  as  such  a  process  enabled  us  to  consider  a  large  number  of 
scatterers  efficiently  and  to  attempt  to  design  signals  and  receivers  which  will 
v/ork  well  on  the  average . 

First,  we  constructed  a  suitable  statistical  model  of  the  random  re¬ 
verberation  return.  Starting  from  physical  considerations,  we  showed  that  the 
return  can  be  reasonably  characterized  as  a  zero-mean  non -stationary  Gaussian 
random  process.  The  correlation  function  of  the  process,  which  is  a  function  of 
the  transmitted  signal  and  distribution  of  the  scatterers  in  range,  is  assumed  to 
be  known. 

Once  this  model  was  obtained,  the  conceptual  path  to  the  optimum  re¬ 
ceiver  was  clear,  and  we  developed  the  equations  which  specify  the  optimum  re¬ 
ceiver  .  The  parameters  in  the  optimum  receiver  depend  on  the  structure  of  the 
reverberation.  The  resulting  receiver  may  be  complex;  frequently,  therefore,  a 
simpler,  sub-optimum  receiver  is  used,  e.g.,  one  which  ignores  the  reverbera¬ 
tion  and  assumes  the  additive  noise  is  white.  In  this  case,  the  receiver  consists 
of  a  filter  whose  impulse  response  is  the  signal  reversed  in  time,  a  detector, 
and  a  threshold  device .  This  type  of  receiver  is  a  conventional  matched  filter 
receiver.  We  have  given  expressions  for  the  performance  of  the  optimum  and 
the  conventional  receiver  for  arbitary  signal  shapes  and  reverberation  charac¬ 
teristics  . 

We  also  considered  the  signal  design  problem .  We  derived  some 
general  properties  that  are  useful  in  choosing  a  suitable  fcrm  for  the  transmit¬ 
ted  signal. 
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To  obtain  a  quantitative  indication  of  the  performance  levels  that  can 
be  achieved  using  various  signal  shapes  and  optimum  or  conventional  receivers, 
we  considered  two  cases  in  detail;  (1)  leverberation  which  is  homogeneous  in 
range  and  results  in  interference  which  is  a  stationary  Gaussian  process,  (2) 
reverberation  which  is  non -homogeneous  in  range  and  results  in  a  non-stationary 
Gaussian  process.  In  both  cases  we  were  concerned  principally  with  a  transmit¬ 
ted  signal  having  a  Gaussian  envelope  and  linear  frequency  modulation. 


B.  CONCLUSIONS  AND  RECOMMENDATIONS 

1.  The  most  effective  way  (within  the  limitations  of  our  model)  to 
combat  reverberation  is  through  proper  .‘?ignal  design.  In  fact,  proper  signal 
design  is  more  important  than  optimum  receiver  design . 

2 .  The  evaluation  of  the  performance  of  the  conventional  receiver 
for  arbitrary  signal  shapes  and  reverberation  characteristics  is  a  straightfor¬ 
ward  calculation.  Finding  the  optimurri  receiver  for  arbitrary  signals,  however, 
appears  both  difficult  and  unrewarding. 

3.  In  many  cases,  the  best  signal  will  not  be  a  Gaussian  pulse.  The 
performance  achieved  using  pseudo-random  waveforms  and  simple  pulse  trains 
should  be  evaluated.  (See  References  14  and  15  for  some  examples.) 

4.  In  some  cases,  the  optimum  receiver  will  offer  enough  improve¬ 
ment  to  warrant  its  complexity.  To  implement  it,  one  needs  to  know  the  statis¬ 
tical  structure  of  the  reverberation  and  thus  needs  some  means  of  measurement. 
Considerable  research  has  been  devoted  to  the  related  measurement  problem  in 
a  radar  astronomy  and  communications  context.  Presumably,  some  of  these 
results  should  be  adaptable  to  the  reverberation  problem . 

5  .  In  most  cases,  the  scattering  function  of  the  reverberation  will 
vary  slowly  with  time.  An  ideal  system  should  have  provision  for  continually 
measuring  the  scattering  function  and  adapting  the  transmitted  signal  shape  and 
receiver  to  the  current  environment. 

6.  We  have  considered  only  detectability  and  time  processing.  The 
problems  of  parameter  measurement  (such  as  range  and  Doppler)  and  the  space - 
time  problem  should  also  be  considered. 
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II.  INTRODUCTION 
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In  the  simplest  detection  problem,  the  signal  returned  from  the  target 
is  completely  known.  In  this  case,  the  receiver  has  available  for  processing  a 
waveform  which  consists  of  either  ambient  noise,  if  no  target  is  present,  or  of 
ambient  noise  plus  the  target  return  signa"  if  a  target  is  present.  The  detection 
problem  consists  solely  of  deciding  which  of  the  two  alternatives  is  correct. 
The  simplest  example  of  this  case  occurs  when  the  ambient  noise  is  a  sample 


function  from  a  white. 


N  9 

Gaussian  random  process  (spectral  height-^  volts'^/cps) 


which  is  independent  from  the  signal.  It  is  well  known ^  that  in  this  case, 
the  optimum  processor  is  either  a  correlation  receiver  or  a  matched  filter 
receiver  and  that  the  performance  depends  only  on  the  ratio  E/N^  (where  E  is 
the  energy  in  the  signal). 


Two  characteristics  of  this  solution  are  of  interest  to  us.  First,  the 
performance  is  completely  independent  of  the  signal  shape.  Any  signal  with  a 
given  amount  of  energy  is  as  good  as  any  other.  Second,  one  can  achieve  any 
desired  performance  level  by  increasing  the  transmitted  signal  energy  to  a  large 
enough  value.  We  shall  see,  however,  that  this  simple  model  does  not  adequately 
describe  the  active  sonai  problem  and  that  these  two  characteristics  cannot  be 
achieved  in  an  actual  situation. 


When  a  signal  is  transmitted  into  the  ocean,  it  encounters  various 
inhomogeneiiies  in  the  medium  and  numerous  objects  which  cause  it  to  be  scat¬ 
tered.  The  return  from  these  various  sources  is  called  reverberation.  (In 
Section  III,  we  shall  construct  a  quantitative  model  of  thi  >  reverberation  return.) 
Since  the  reverberation  return  is  caused  by  the  signal,  it  is  clear  that  the  statis¬ 
tical  characteristics  of  the  reverberation  "noise"  are  not  independent  of  the  sig¬ 
nal  shape.  One  would  suspect,  therefore,  that  the  receiver  performance  will  no 
longer  be  independent  of  the  signal  shape.  Moreover,  increasing  the  transmitted 
signal  power  will  increase  the  level  of  the  reverberation  return.  Thus,  increas¬ 
ing  the  transmitted  energy  may  be  an  inefficient  way  to  combat  reverberation. 

The  problem  bere  is  one  of  target  detection  in  the  presence  of  inter¬ 
ference  which  depends  on  the  transmitted  signal.  The  basic  ideas  involved  in 
our  solution  are  reasonably  straigiitforward,  but  the  manipulations  necessary  to 
obtain  a  quantitative  solution  are  some"''hat  involved.  To  illustrate  some  of  the 
basic  concepts  involved,  we  will  consider  a  simple  example  of  target  detection 
in  the  presence  of  interference.  The  model  does  not  represent  a  realistic 
sonar  problem,  but  is  only  a  tutorial  exampi<=*. 

/We  assume  that  the  reader  is  familiar  with  the  application  of  statistical  detec¬ 
tion  theory  to  receivers  operating  in  an  additive  Gaussian  noise  environment. 
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Consider  the  simple  problem  shown  in  Figure  i.  The  signal  returned 
from  the  target  is  Sjj(t).  The  signal  returned  from  the  interfering  object  is 
Si(t).  In  addition,  an  additive,  ambient  noise  n(t)  is  present.  The  received 
waveform  is  a  sum  of  these  three  terms. 


r(t)  =  S^(t)  +  Sj(t)  +  n(t)  (ll-l) 

As  a  special  case,  let  us  assume  that  Sjft)  is  identical  to  S^ft)  except 
for  a  time  delay  and  an  attenuation. 

Thus, 


Sj(t)  =  aS^(r  -  T) 


(11-2) 


where  t  is  assumed  known.  If,  in  addition,  the  value  of  "a"  is  known,  the  solu¬ 
tion  is  simple.  One  subtracts  out  the  interfering  signal  and  then  uses  the  usual 
matched  filter  receiver.  The  receiver  structure  for  this  simple  case  is  shown 
in  Figure  2.  To  make  the  problem  more  realistic,  we  must  include  some  uncer¬ 
tainty  to  the  interfering  signal.  Therefore,  we  assume  that  the  attenuation  is  a 
Gaussian  random  variable  with  zero-mean  and  variance,  Oa  . 


We  now  have  a  familiar  two-hypothesis  problem. 

Under  Hq,  no  signal  present,  the  mean  of  r(t)  is  zero  and  the  covar¬ 
iance  function  is: 


N 

R(t,u)  =  E  rr(t)r(u)]  =  6  (t-u)  +  a|  S^(t-T  )S^(u-t  )  (II-3) 

where  we  assume  the  additive  noise  is  a  sample  function  from  a  Gaussian  process. 

Under  H^,  signal  present,  the  mean  of  r(t)  is  and  the  covariance 

function  is  unchanged. 


/  This  example  is  based  on  a  similar  example  in  an  unpublished  memorandum 
byW.  M,  Sieben. 


4 


Arthur  Q.Hittle.Jnr. 


FIGURE  2  OPTIMUM  RECEIVER- INTERFERING  SIGNAL 
KNOWN  COMPLETELY 
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The  optimum  detector  consists  of  a  correlation  operation  on  the 
received  signal  r(t) 


.r„ 


X  =  1  r(t)  q(t)dt 
0 


(II-4) 


where  q(t)  is  the  solution  to  the  integral  equation: 


q(t)  =  J  Q(t.u)S^(u)  du  0<t<T 
0 


(11-5) 


and  the  function  Q(t,  u)  is  the  "inverse  kernel"  and  satisfied  the  equation: 


6(t-z)  =  j  Q(t,  u)  R(u,  z)du 
0 


(II -6) 


This  solution  is  just  a  special  case  of  the  nonwhite  noise  problem. 
One  can  verify  by  direct  substitution  into  Equation  II -6  that 


Q(t.U)  = 


6  (t-u)  - 


^a 


N  N 


o 


o 

2 


+ 


S^(t-T'S^(U-T)  (II-7) 


Then  substituting  Equation  11-7  into  Equation  II-5,  we  obtain: 


q(t)  = 


N- 

o 


N 


a  -  E  -f  N 


j  S^(^)Sjj(u-')du 


N 


a  o 


(II -8) 


6 


Arthur  2).l.tub.9nr. 


where 


j  S  (u)S  (u-T  )du 

v,,<n-^-4 -  (n-9) 

represents  a  normalized  correlation  between  the  desired  signal  and  the  interfer¬ 
ing  signal,  Clearly,  O^Y  (t)^1. 


The  optimum  receiver  consists  of  two  parts,  as  shown  in  Figure  3. 
One  part  is  the  usual  correlation  operation;  the  second  pait  is  a  partial  subtrac¬ 
tion  of  the  interfering  signal. 


S(jIt)  2/No 


FIGURE  3  OPTIMUM  RECEIVER:  INTERFERING  TARGET 
WITH  RANDOM  AMPLITUDE 
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The  optimum  receiver  has  a  very  good  intuitive  interpretation.  In 
Figure  2,  we  saw  r'nat  if  the  interfering  signal  were  known  exactiy,  we  couid 
subtract  it  out.  A  logical  approach  in  the  presence  of  some  uncertainty  might 
be  to  look  at  the  received  waveform  r(t)  and  estimate  what  the  interfering  signal 
is.  We  could  then  subtract  this  estimate  and  pass  the  result  into  the  normal 
white  noise  detector.  We  will  now  demonstrate  that  this  logical  approach  is 
exactly  what  the  optimum  receiver  is  doing. 

The  only  unknown  quantity  in  the  interfering  signal  is  the  amplitude 
"a”.  One  can  show  that  if  the  received  waveform  is 

r(t)  =  aS^(t-T)  +  n(t),  (II-IO) 

then  the  most  probable  value  of  "a"  and  the  minimum  variance  estimate  of  "a" 
are  identical .  The  estimate,  a,  is  given  by 

o^  X 

0 

2 

The  equation  describing  the  receiver  shown  in  Figure  4  is: 

X  =  J  [  r(t)  -  a  S^(t-T )]  S^(t)  dt  (IM2) 

0 


Substituting  Equation  11-11  into  Equation  11-12,  we  obtain: 

T  CF®  Y  T 

X  =  j  ^  -  ’  J  r(u)S^(u-T)du  (11-13) 

OO  "^-I+o^eO 

2  a  ^ 

We  see  that  the  two  receivers  are  equivalent.  Thus,  the  optimum 
receiver  does  exactly  what  one  might  expect. 


A  complete  measure  of  performance  is  the  ratio  of  the  square  of  the 
mean  of  X  under  hypothesis  H ^  to  the  variance  of  X. 


Thus. 


d"  H 
o 


r  '»  3 

\E  [X  :  Hi  ]  ; 

Var  X 


(11-14) 
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One  can  show  that 


o 


1  - 


i£fc 


a. 


1  + 


2E 

N 


(IMS) 


We  observe  that  the  first  term  is  simply  the  usual  white  noise  result. 
The  second  term  represents  the  degradation  due  to  the  interfering  signal.  The 
magnitude  of  this  degradation  depends  on: 

1.  (T)  :  the  correlation  (or  similarity)  between  the  desired 

signal  and  the  interfering  signal . 

the  strength  of  the  interfering  signal 

the  energy-to-noise-density  ratio. 

If  any  of  these  are  small,  the  effect  of  the  interfering  signal  will  be  small.  The 
signal  design  problem  in  this  case  is  simply  choosing  a  signal  shape  so  that 
Y  ,.  (T )  is  small.  For  fixed  t  and  no  constraint  on  peak  transmitter  power, 
tie  solution  is  simple. 

Let 


s,(t) 


271" 

T 


0<  t  ^  J 


0  Elsewhere 


(11-16) 


In  general,  for  the  cases  of  interest,  the  solution  will  be  more  com¬ 
plicated  . 


Now,  let  us  assume  that  the  designer  is  unaware  of  or  chooses  to 
ignore  the  interfering  target.  He  would  then  use  a  conventional  matched  filter 
which  is  now  nonoptimum.  One  can  verify  easily  that 


cf 


2E 

N 

o 


1 

X 

1  +  Y,.^T) 


(11-17) 
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!t  is  convenient  to  compute  a  cjuantity,  '  the  deg^adsfion  due  to  inter¬ 
ference.  "  Taking  the  log  of  Equation  11-15,  we  have: 


log  =  log 

o 


r  o» 

N  a 

+  log  1-V=^,(T)  - -° 

1  +  ^  a 
M  a 

^  o 


(11-18) 


Since  the  first  term  is  caused  by  white  noise,  the  degradation  due  to 
interference  is  just  the  magnitude  of  the  second  term. 


Similarly,  from  Equation  11-17, 


logd^  =  log  ^  -  log  1  +  Y 


»  (T)  .  ^  0=  1 

di  '  '  N  a 

o  J 


(11-19) 


The  degradation  for  the  conventional  and  optimum  receiver  cases  is 
shown  in  Figure  5. 

We  observe  that  in  both  the  good  performance  region  (Y*  -*0)  and  the 
bad  performance  region  (y^  -*  1) ,  the  optimum  receiver  is  not  much  better  than 
the  conventional  receiver.  We  can  see  that  this  result  is  intuitively  logical  by 
looking  at  Equation  II -8. 

As  Y*  -•  0,  the  coefficient  of  the  second  term  approaches  zero. 
Physically,  Y*  =  0  means  that  the  desired  signal  and  the  interfering  signal 
are  orthogonal.  Thus,  the  interfering  signal  causes  no  outpu:  in  the  correlation 
detector.  Clearly,  a  signal  that  causes  no  output  cannot  affect  the  performance, 
and  there  is  no  reason  to  modify  the  detector. 

At  the  other  extreme,  as  Y^  -  1,  the  modifying  term,  S{j(t-T ),  looks 
more  and  more  like  the  original  term.  In  the  limit,  y*  =  1,  S^(t-T)  =  S^(t) 
and  no  modification  is  necessary. 

This  simple  example  illustrates  many  of  the  important  features  of  the 
actual  reverberation  problem.  W'e  may  summarize  these  briefly: 


1 .  The  optimum  detector  tries  to  subtract  out  the  interfering 
signal.  Since  it  does  not  know  the  signal,  it  uses  the 
received  waveform  to  estimate  the  interfering  signal  and 
then  subtracts  out  this  estimate. 


2.  The  choice  of  sigiiais  affects  tiie  performance  appreciably. 
The  goal  of  the  signal  designer  is  simple  to  state:  make  the 
signal  returned  from  the  desired  target  orthogonal  to  all  of 
the  interference.  We  will  see  that  in  most  realistic  situa¬ 
tions,  as  we  would  expect,  this  is  difficult  to  do. 

3.  To  design  an  optimum  detector,  requires  some  knowledge 
of  the  interfering  signal.  In  this  simple  case,  we  knew  the 
shape  and  the  probability  density  of  the  amplitude.  We 
would  expect  that  as  knowledge  of  the  interference  decreased, 
the  improvement  gained  would  decrease.  In  many  cases, 
because  of  a  lack  of  knowledge  or  in  order  to  simplify  the 
resulting  equipment,  one  uses  a  "conventional"  receiver. 

We  observed  that  for  certain  parameter  ranges  there  was 
not  too  much  difference  between  the  conventional  and  the 
optimum  receiver.  We  will  find  f'lat  in  many  cases  of 
interest  proper  signal  design  is  much  more  important 

than  the  difference  between  a  conventional  and  an  optimum 
receiver. 


B.  PROBLEM  FORMULATION 


There  are  important  differences  between  the  simple  problem  discussed 
above  and  the  reverberation  problem: 


1 .  Instead  of  one  interfering  target,  there  is  a  large  number 
of  interfering  returns  from  reflecting  objects. 


2.  The  returned  signal  from  the  target  is  a  band-pass 
waveform.  It  has  a  random  phase  angle  which  must  be 
taken  into  account. 


We  shall  see  that  these  differences  take  us  from  a  tutorial  exercise 
to  a  reasonably  good  model  of  an  active  sonar  in  a  reverberation  environment. 
The  cost  of  this  transition  is  a  great  increase  in  the  complexity  of  the  calcula¬ 
tions.  It  is  important  to  emphasize  that  the  concepts  in  an  actual  sonar  problem 
are  identical  to  those  in  the  preceding  example. 


Our  model  of  the  reverberation  problem  is  shown  in  Figure  6. 
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SCATTERERS 


TARGET  ( DELAY  Tq  ; 


DOPPLER  L)o) 


SIGNAL  RETURN:  A  Re  [Tlt-Tp 

+  REVERBERATION  RETURN:  NrH) 


RECEIVER 


AMBIENT  NOISE;  Na  (t) 


FIGURE  6  CHANNEL  MODEL 


The  tiansmitted  signal  is: 


S^(t) 


=  Re 


r  1  r 

f(t)e  =  Re  u 


t  ■ 


u(t)  COS(UL’^t  + 


0(t)  (n-20) 


where:  , 

f(t)  is  the  complex  envelope 

u(t)  is  the  actual  envelope 

0(t)  is  the  phase 

/  For  a  discussion  of  the  complex  representation,  see  Woodward  (Reference  1) 
and  Hel  Strom  (Reference  2). 


14 


Arthur  Zl.lLittle.inr. 


/ 


I  The  returned  signal  R(t)  =  Re 


1 


% 


n  (t)e 


mt\  t 


consists  of  the  following  three  parts: 


I .  The  reverberation  return  due  to  the  collection  of  nonstationary 
scatterers: 


Nj^(t)  =  Re 


r  ' 

n  (t)e 
r 


(11-21) 


2.  The  return  due  to  additive  Gaussian  noise; 


N^(t)  =  Re!  n^(t)e 


r  '"c'l 


(11-22) 


I 

I 

J 

i 

I 

i 

I 


I 


3.  The  return  due  to  a  target.  This  is  an  attenuated  and  phase- 
shifted  replica  of  the  transmitted  signal . 

[+  jS  +  ji)  t  +  7U  t  I 

f(t)-T^)e  J  (11-23) 

where  B  is  uniform  [  O,  2tt]  and  unknown.  The  first  two  parts  of  the  return 
are  always  present.  The  third  part  is  only  present  if  a  target  is  present. 

As  before,  we  are  concerned  with  deciding  whether  or  not  a  target  is 
present.  In  other  words,  we  want  to  decide  between  one  of  two  hypotheses: 


Hq  (no  target): 


r  (t)  =  n  (t)  +  n  (t) 
ora 


(target  present): 


ju)  t  +  i6 

r (t)  =  f  (t  -  Tj^)e  +  n^(t)  +  n^(t) 


If  n  (t)  is  a  complex  Gaussian  process,  the  solution  to  the  problem 
can  be  ejqjressed  in  terms  of  an  integral  equation.  Assuming  the  integral  equa¬ 
tion  can  be  solved,  one  can  then  construct  the  optimum  receiver. 
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7^ 


III.  MODEL  FOR  NON-UNIFORM  SCATTERERS 


In  this  section,  we  develop  a  model  for  the  reverl>eration  return. 

Our  approach  is  a  genet alization  of  that  in  Reference  3  (more  readily 
available  References  are  4  and  S'!. 

The  transmitted  signal  is: 

r  1 

S^(t)  =  ReLf(t)e  (in-l) 

We  are  concerned  with  f  (t)  the  complex  envelope. 

?(t)  =  (III-2) 

'Fhe  complex  envelope  returned  from  an  individual  scatterer  (the  n 
scatterer)  is: 

jX'  f 

S  (t)  =  Z  e  T(t-t)  (111-3) 

n  n  n 


where: 

Z  is  a  complex  number  which  is  the  magnitude  and  phase  of 
the  echo  (i.e.,  'strength'*  of  echo), 

th 

x  is  the  Doppler  shift  due  to  the  radial  volccity  of  the  n 
scatterer,  and 

t^  is  the  delay  due  to  die  position  of  the  n  scatterer . 

The  entire  complex  envelope  due  to  reverberation  is: 


Etx  t 

n  ^ 

Z  e  "  f  (t  -  t  ) 
n  n 

all  scatterers 

We  make  the  following  assumptions  regarding  the  scatterers: 


(1II-4) 
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Assumption  1 

The  distribution  along  the  path  obeys  a  non -homogeneous  Poisson  law. 
(See  Reference  6.) 

The  probability  that  a  scatterer  exists  in  the  time  interval.  It  is: 


Pr  1  i  event,  t  <t<t  l=a(t  )dt 
L  a  aj  ' 


(III-5) 


All  of  the  properties  of  a  stationary  or  homogeneous  process  can  be 
extended  easily  to  a  non-homogeneous  process.  The  two  properties  that  we  will 


use  are: 


aT  \n  T 

a  (x)dx  I  exp  -  T  a  {x)dx 

_  T  /  J-T 


Property  I:  Pr(n  events  in  interval  [-T,T])  = 


Property  II:  Given  that  n  events  occur  in  the  interval  [-T.T]  ,  the  joint 
probability  density  of  their  occurrence  times  is  given  by  the  expression 


(111-6) 


I  ,  .  ,  Si'  ^2 . ^n] 

I  n  events  in  (-T,T) 


r  T  1^ 

Ja(x)dx 
T 


a(t,)a(t_) . a(t  )  (HI-?) 

l  z  n 


Assumption  2 

The  velocity  of  each  scatterer  is  a  random  variable.  Velocities  of 
different  scatterers  are  indepiendent  random  variables . 


The  probability  density  governing  the  scatterer  velocity  is  time 
dependent.  The  probability  that  a  scatterer  occurs  in  the  intervalj^t  ,t  -tdt  J 
and  has  a  velocity  (frequency  shift)in  the  range  can^ro  ^ 

described  by  a  joint  density. 


r  1  1, 


p  (t  ,  Ju  )  =  p  I  (o)  1  t  )  p  (t  )  =  p  ,  )  '  --IW-  - -  iIII-S) 

^t,  ju'a  r  t  l>  a' ^t  a'  ^aj't'  1  a  r^a{y>dx 
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AssumDtion  3 


The  strength  nf  each  scBtterer  is  a  random  variable.  Strengths  of 
different  scatterers  are  independent  random  variables.  The  strength  is  independ¬ 
ent  of  both  position  and  velocity.  We  denote  this  probability  densiry  by; 

Pz  =  Pz  (2)  =  . . .  p^  (Z)  (III-9) 

1  2  n 

Using  these  three  assumptions,  we  want  to  find  the  correlation  function 
of  the  reverberation  return.  For  simplicity,  we  will  assume  E  [Zjjj  = 

The  complex  envelope  of  the  returned  reverberation  signal  is; 

m  =  0 

Then,  the  correlation  function  is; 

R  tg  )  =  i  <  r  (t^)  r^tg)  >  (III- 11) 


An  easy  way  to  find  this  is  to  assume  there  were  n  scatterers  in  the 

interval  (-T,T^  If  we  denote  the  conditional  correlation  function  based  on  this 

assumption  by  R  (t  .  t^),  then: 
nos 


R  ^n^^o’  ^8^  ‘  events  in  ('T,T)] 

n  =  0 

Using  Equation  III-6,  this  reduces  to; 


n  =  0 


RJt 


n  a 


~  rT  le  /.y 

I  a(x)dx  exp  -  /  a(x)dx 

U-T  J  J-T 

n ! 


(Ill- 12) 


(III- 13) 


/  We  use  the  symbol  E  to  denote  the  expectation  of  a  rand"'m  variable. 
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i 


I 


1 


g 

Jist 


Now  we  find  R  ft  .  t .  V 
rr  a  »' 


using  equation  m-/  tnrough  iii-ii,  we  may  wrilc; 

T  T 


dt, . .  .dt  a(t,). .  .a(t  ) 
1  n  1  n 


/•/  t, 

•'•'11  n 


■  30  -30 


t  (u)  :  t  )dt«, . .  .duu 
n  n  n  1  n 


.OS  -00 


n  n 


e 


i  =  1  k  =  1 

Several  observations  simplify  this  expression: 

or  OD 

1.  / /p^  (Z,)P  (Z^'Z,ZJ  =E|Zj|'  l  =  k 

0  ifk 

This  reduces  the  double  summation  to  a  single  summation. 

2.  For  i  =  k,  the  exponential  terms  reduce  to: 


+i*,(ta-tg)  +rx..T 

e  =  e 


(where  T  =  t  -t.) 

a  B 


(III- 14) 


(in -15) 


(in- 16) 


/  This  approach  is  similar  to  Sec,  7-4  in  Reference  7.  This  section  follows  the 
original  work  of  Rice  (Reference  S). 
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Now  consider  the  term  in  the  series  for  i  =  q.  After  performing  the 
integration  with  respect  to  Z_,  the  multiple  integral  becomes; 


T  T 


'  f  f 

— Tf - —  /  . . . /  dt, . . .dt  , 

ra(x)dx  "  J  J  ^ 

J.T  -T  -T 


dt  . .  .dt  a(t  ). .  .a(t  )a(t  ). .  .a(t  ) 
q+1  n  1  q-1  q+1  n 


q+11  q4- 


...p  ,  (tn  ;t  )  di’, . .  .dju  dju  , . .  .djj 

1 1  n  n  1  q-1  q+1  n 


(Zi)...Pz  (z,.i>Pz 

d  {  I  q-1  ^  q+1  n  q-1  ‘ 


•  00  -00 


f  a(t  )dt  fp  (a'  It  )daj  •  E  (\Z  -t  )r(t.-t  )e  ^  (III-17) 

J  'q'  q  J  ql  q'  q  |!  q!  f  '  a  q'  '  6  q 


The  integrals  with  respect  to  variables  other  than  t^  and  ju^  are  straightforward . 

We  observe  that  the  integral  with  respect  to  x  is  in  the  form  of  a  condi¬ 
tional  characteristic  function: 


M,  (-t  ) 

X'  •  t  q 

q  q 


/+j'XaT 

p  1  (x  1 1  )e  *  dx 
X  t  '  ql  q-  c 


(III- 18) 


Assuming  that  T  is  large,  we  can  neglect  end  effects. 

Since  each  term  in  the  sum  is  the  same,  we  have: 


rx  ,  .X  i  n  E  i  z!“ 

R  (t  .t  )  =  t  - 

n  a  3  fi 

I  a (x)  dx 

•/ X 


-J —  f  a(t.)M  ..  )  (III-19) 

j  4  .  .  4  -  q  -  q 

-T 
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'cr.ote  the  integral  by  !  (t  .t^).  Suhstitutirg  Rquarion  111-15  into  Equation  III-9,  we 


have: 


n  =  0 


a(x)dx  exp  -f  a(x)dx  „  gj-,  ^  2-| .  ,(,^  (5) 

n !  '  /-T  " 

J  a(x)clx 


(III-20) 


This  reduces  to: 


JU 

T 

R(t^.tg)  =  iE[.Z|-]l(t^.tg)exp-  Xr 


-]n-l 

a(x)dx 

ll _ J_ 

(n  -  1) ! 


(111-21) 


We  observe  that  the  sum  cancels  the  exponential  term  preceding  it, 
Equation  ?n-21  becomes; 


R(t^.tB)  =  iE[!z!^]  I(t^.tg) 


(III-22) 


where 


‘  K’  's>  '  X®  :x“a''e:x>  ‘  <‘o'  '  <‘s  '  ^ 


(III-23) 


Equation  III-23  is  valid  for  the  case  defined  by  the  original  assumptions, 
and  it  can  be  written  in  several  different  ways .  The  two-dimensional  correlation 
function  of  a  signal  is  defined  to  be'X 


/~  T  T 

(t+2) 


-jiJt . 
e  dt 


an-24) 


.V.  T  T  1  /* 

f  (t-2)  f  <«+2>=  2;;  J 


jj)  e  do; 


(ni-25) 


By  letting 


t  +  t^ 

t  =  — - - x,  T  =  t  -  t  ,  Equation  111-25  becomes; 


/  See,  e.g. ,  Woodw'ard  (Reference  1)  or  Siebert  (Reference  9). 
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f  (t^-x)f  (tj-x)=—  (tg-V»=)e 


* 


\ 


du) 


(III-26) 


Substituting  Equation  !TT-26  into  Equation  III-22,  we  have; 


J |*a(x)  |,(>„-tg:x)  ej(tg-t^;a,)  exp +ja. - xj 


2 

(III-27) 


/■ 


But  /  a  (x)  (t^  -  tgtx)  e  •’^dx  =  S  (  -oj;  t^  -  tg) 

q' 


(III-28) 


where 

S(Vi,V2)=//e  ^,(x,,X2)dx,dX2 


[ 


die  characteristic  fctn 
the  joint  density 


.  of  j 


(in -29) 


2  ^ 

‘3>  '  ^  V  ‘s’  ^  <‘s  ■  ‘a  “P 

(III -30) 

The  total  interfering  noise  is  the  sum  of  the  reverberation  noise  and  die 
ambient  noise. 

The  total  correlation  function  is:  ^ 


=  +  (t„.U 


r  a 


o  'a  B'  n  'a’  B' 
r 


(III-31) 


For  a  large  number  of  scatterers,  one  can  show  that  the  interfering 
noise  approaches  a  non -stationary  Gaussian  process. 


/  We  assume  the  ambient  noise  is  a  real  white  Gaussian  process  with  double - 
sided  spectral  height  No/2 . 

//  The  technique  is  similar  to  that  of  Reference  4. 
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I 


I 

I 

I 

I 

I 


The  simplest  case  is  when  the  distribution  in  range  of  the  scatterers 
is  uniform  and  the  velocity  density  is  uniform. 

Let  the 

a(x)  =  Y  -T<x<T  (III-32) 

Physically,  v  is  the  number  of  scatterers  per  unit  interval  (time).  Now  let 
T-=o  . 

Then,  from  Equations  III-22  and  III-23 

+  00 

R(t  .tJ  =  t  yE[Iz!^]M,,  (t  -U  fl  (t  -X)  r (t„-x)dx  (III-33) 

O  C  'Jt)  OL  D  f  CL  a 

q  J 


or 


i 


T 

! 


I 


I 

I 

I 

I 

I 

I 

I 


letting  u  =  ■  X 


a  D 


we  have: 


+  a> 


,2.. 


R(T) ye[;  zriM^  (T)  pill)  r (u+T)du 

q  J 


(in-34) 


But  the  integral  is  just  Rj(t) 


R,(T) 


f  (u)  (u  +  t)  du 


(in-35) 


Then, 


(T)  =  ^  Y  E  [I  Z!^]  (T)  R^(t)  h  ^  (T)  R^(t)  (III-36) 

r  '"q  '"q 

Our  principal  results  in  this  section  are  Equation  in-30,  the  correlation 
function  for  the  reverberation  return  for  the  case  of  non-uniform  distribution  in 
range  of  the  scatterers,  and  Equation  111-36,  the  correlation  function  for  the  rever¬ 
beration  return  for  the  case  of  uniform  distribution  in  range  of  the  scatterers . 

Since  the  process  is  Gaussian  it  is  completely  characterized  by  its  correlation 
function.  In  the  next  section,  we  use  this  to  develop  the  optimum  receiver  structure. 
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IV. 


DERIVA'nON  OF  OPTIMUM  RECEIVER  cTpiirrTipp 


4^1  t 


In  this  section,  we  derive  some  general  results  regarding  receiver 
structure  and  performance  that  we  will  need  for  our  specific  problem. 


A.  STRUCTURE/ 

We  first  derive  the  structure  of  the  optimum  receiver  and  the  conven¬ 
tional  receiver.  Then,  we  find  expressions  for  their  performance. 

As  pointed  out  in  the  introduction,  the  optimum  receiver  solves  the 
hypothesis  testing  problem.  One  can  demonstrate  that  for  many  interesting 
criteria  (e.g. ,  Bayes,  Neyman  -  Pearson,  Minimax)  the  solution  reduces  to  one 
of  forming  the  likelihood  ratio  and  comparing  it  with  a  threshold./  The  value  of 
the  threshold  will  depend  on  the  decision  criterion  and  relative  costs.  We  will 
not  concern  ourselves  with  choosing  a  specific  value  of  the  threshold  but  only 
with  forming  the  test  statistic. 

We  are  concerned  with  detecting  a  signal  which  is  a  member  of  the 
ensemble  S(t, S),  where: 


S(t,9)=  Re  [s^(t)ej^e^  J  (IV-1) 

where  ~  ^  0  <  3  <  ’rr  (IV-2) 

and  S  (t)  =  Af(t-Tj)  e  (IV- 3) 

d 

To  form  the  likelihood  ratio,  we  expand  the  complex  envelope  r  (t) 
using  a  Karhunen- Loeve  expansion.// 


/  Helstrom  (Reference  2),  C3iapter  3. 
//  See  Reference  7,  Chapter  6. 
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The  orthonormal  functions  of  this  expansion  are  the  eigenfunctions  of 
the  integral  equation: 


=  /  Rj.  <ta«t3)  «>n(t3) 


dt; 


-T 


where  Kj.  (tg^,  ta)  is  the  complex  covariance  of  rQ(t). 


(IV-4) 


=  i  E  [  (,3)] 


(IV- 5) 


Tlien, 


r(t)  =  E  rn0n(t) 


T 

=  J  r  (t)  0;  (t)  dt  =  Xn  +  j 


Vn 


(IV-6) 

av-7) 


We  now  want  to  find  the  statistics  of  and  under  the  two  hypotlieses. 

Under  the  no-target  condition  (hypothesis  Hq),  the  complex  envelope  of 
the  received  signal  is 


O  X  cL 


(IV-8) 


Under  the  target  present  condition  (hypothesis  H^),  it  is 


T^it)  =  S^(t)e^^  +  n^(t)  +  n^(t) 


(IV-9) 


Clearly,  x^,  and  y^^  are  Gaussian  random  variables  under  either  hypothesis. 
Under  hypothesis  Hq, 

®  ["n]  = 
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while  under  hypothesis  Hj, 


J 

L  -T  “* 


(IV- 11) 


where 


S  = 
n 


(t)  (t)  dt 


(IV- 12) 


The  covariances  are  independent  of  the  hypothesis  and  are: 


E  =  f  d-J  f 

LnnmmJ  j  J 


dv  E  f  (u)  r  *  (v)  1  0*  (u)  0  (v) 

L  J  n  m 


-T  -T 


/•T  /’T 

I  du  I  dv  R  (u,v)0*(u)0  (v) 

II  r  n  m 

J  J  o 


-T  -T 


Using  Equation  IV- 4,  we  have: 


(IV- 15) 


efc  -r  )(r*  -r*  )1  =  2\  f  du  0* (u)  0  (u)  =  2X  6 

LnnmmJ  ^  J  ^  ^  mi 


(IV- 16) 


Similarly,  using  the  propeities  of  the  expansion,  we  have: 


E  [  1-  r  1  =  0 
L  n  mj 


Since  r  =  x  +  j  y 
n  n  n 


(TV- 17) 


(IV- 18) 
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i 


r  1  r  _ 

we  have  E  (x  -  x  )  (x  -  x*  )  =  *.  ^^  ;  E  (y  -  v  )  (y  -  y*  )  I  =  ^ 

1_  n  n  m  in  J  n  nm  ^  n  'n  m  m  J  n  nm 


(IV- 19) 


[< 


and  E  (x  -  x  )  (y  -  y  ) 

L  n  n  m  m 


(IV-20) 


Therefore,  the  x^  and  yj^  are  statistically  independent  random  vari¬ 
ables  .  Now  write  the  likelihood  ratio: 


P(r  I  Hj) 

W\%) 


(iv-2n 


Now 


'(?l  Hj)=j 


2tt 


p(S)d8 


K 

TT 

i=J 


2ttXi 


K 

e*P  -Z 

n=l 


(x  -E(x  ))^  +  (y  -E(y 
n  n  -^n  n 

2X 

n 


(IV-22) 


or 


,r  IH^) 


r 


p(S)d0l 


exp 


2X 

n 


(IV- 23) 


and 


(IV- 24) 


We  have: 


1  ^ 

2Rej 

r  *  -j0 

2 

-  s 

1  -  A<VT>  ^  i 

y  ' 

L  n  n  . 

n 

J  2r,  ^ 

2X 

0  j 

n=l 

n 

av-2s) 
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(IV- 26) 


A  =  T  "  ^ 

n=l  n 


where  is  real  and  positive. 
Then,  we  nrsav  write: 

■  ^2”  1 
“  =  2?  J  <19  exp  -  j 

o 


cos  (2l  -  3)  ^ 


rv  ,  ^ 

1  n 

n=i  n 


(IV-27) 


But  the  e'Dression  in  tlie  bracket  is  just  iQ(Aj^)  where  is  a  modified  Bessel 
function  of  the  first  kind  and  order  zero.  That  is 


‘o  T. 


IM 

2X 


Now  Iq(*)  is  monotone  for  positive  and  negative  arguments  and  :^^ym- 
nietric  around  the  origin.  Letting  K  -  ®,  we  see  that  an  adequate  statistic  is: 


K 

* 

s  v 

fT 

L 

n  n 

\ 

= 

/  q*  (t)  r  (t)  dt 

n-l 

n 

-T 

(IV-28) 


where  q*  (t)  is  the  solution  to  the  integral  equation. 


s,  :t ) 

d  a 


J  o 


a’V  *3 


(lV-29) 


The  second  term  in  Equation  IV-27  is  incorporated  into  the  threshold. 
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TTie  desired  ootinium  operarjo.os  ore  shown  in  Figure  7.  7 he  desired 
operations  can  be  realized  pirysically  by  passing  r(t)  through  a  narrow-band 
filter  whose  complex  impulse  response  is 


h  (T)  =  q*(T-T)  -T<T<7 

opt 

h  (T)  =  0  elsewhere  (IV- 30) 


and  detecting  the  output  envelope.  This  realization  is  shown  in  Figure  S. 

We  now  want  to  consider  three  cases: 

1.  The  scatterers  are  uniform  in  range.  This  leads  to  a 
stationary  process  for  n-pCf)  (see  Equation  IV-36).  One 
solves  Equation  IV-29  to  find  the  optimum  receiver. 

2.  The  scatterers  are  non-uniform  in  range.  This  leads  to 
a  non- stationary  process  for  ny(t).  One  solves  Equa¬ 
tion  IV-29  to  find  the  optimum  receiver. 

3.  The  reverberation  return  is  ignored  in  finding  the  opti¬ 
mum  receiver.  We  assume: 


=  if 

o 


(IV- 32) 


Fhis  is  called  a  conventional  receiver.  We  then  investi¬ 
gate  the  performance  of  the  conventional  receiver  for  non- 
uniform  and  uniform  scatterer  distributions. 


4  In  the  next  section,  we  consider  cases  (1)  and  (3). 


q*(t) 


FIGURE  7  DESIRED  OPERATIONS  (COMPLEX  ENVELOPES) 


NARROW- SAND 
FILTER 


FIGURE  8  PHYSiCAi.  REAUZAT’^  N 
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B.  PERFORMANCE  OF  OPTIMUM  RECEIVER/ 


Our  decision  is  based  on  comparing  the  magnitude  of  A  with  some 
threshold  ''>0- 

We  must  consider  the  statistics  of  Ia  j  under  hypothesis  Hg  and 

hypothesis 


Pr  [False  Alarm]  =  J' 


p  (A)  dA  =  P 

O  r 


(IV-33) 


X  +  jy 


(IV- 34) 


X  -  Re 


(t)  rvt)  dt 


aV-35) 


Var  X  =  Var 


r  - 

•  y  =  J  q 


(t>  Sj(t)  dt  =  d 

d  o 


aV-37) 


P,(A)  =^e 
o 


(iV-38) 


/  Our  derivation  is  rather  sketchy.  The  details  are  on  pp.  149-156  of  Helstrom 
(Reference  2). 
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where 


V  =  Im 


(t)  dt 


(IV-36) 


Hu. 


Substituting  Hountion  IV-38  into  Hcjustlon  IV~33  snd  svElusting,  wc?  hsvs! 


P 


F 


=  exp 


Similarly,  one 


can  show,  under  h^jpothesis  H^,  that 


{IV-39) 


2 

X  =  d  cos  0 
o 


y 


d  sin  S 
o 


(IV- 40) 

(IV-41) 


and  the  probability  of  detection  is; 


A  . 

Pr  [detection]  =  )  =  Pq  (IV-42) 

o 

where  Q(a,b)  is  Marcum's  Q  function.  (See  References  8,  11,  and  12.) 


Q(a,b) 


OD 


X  exp  - 


b 


/  2  2, 
(x  +a  ) 


I  (ax)  dx 
o 


ra3serve  that  from  Equations  IV-39  and  IV-42,  we  may  write: 


P_  =  Q(d,-21nP_)  (IV- 42a) 

U  O  r 

2 

One  can  plot  Pp  vs.  d^  as  a  function  of  Pp.  This  curv^e  is  shown  in  Figure  9./ 

When  the  correlation  function  is  stationary,  there  is  a  sim.ple  expres¬ 
sion  for  d^  in  terms  of  the  various  spectra.  For  simplicity,  assume  that  the 
observation  interval  is  infinite. 


/  Our  figure  is  similar  to  Figure  V-2,  p.  155,  Reference  2.  An  earlier  refer¬ 
ence  is  Reference  10. 
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FIGURE  9  RECEIVER  OPERATING  CHAR^iCTERISTIC: 
OPTIMUM  RECEIVER 
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Tlien,  Equation  IV- 29  becomes; 


Ti 

/ 


^  I  R(t  -t-)  q  (t^)  dt- 

d  a  I  a  ;  : 


(IV- 43) 


and  Equation  IV- 37  reduces  to; 


CD 

■o  '/ 


(t  )  S.(t  )  dt 
a  d  a  a 


(IV-44) 


Equation  IV- 43  can  be  solved  using  Fourier  transforms.  Transforming, 


we  have; 


Q(0  = 


S^(x) 


S,<x) 

S  (1.)  ‘  N  +S  (a) 

I\  )  u 

r 


(IV- 45) 


and  using  Parseval's  theorem: 


2  r  I 

°  '  J  V^> 


d'-t 

2^ 


(IV-46) 


In  the  absence  of  reverberation,  the  receiver  is  specified  by  Equation 


IV- 32  and 


o  N  J 

o  •' 


2E 

S,(t)S‘(t)dt  =  ^ 


(IV- 47) 


Observe  that  d  is  just  the  S/N  ratio  at  the  receiver  output, 
o 


* 

- 1 


w 

^1 
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C.  PERFORMAN’CE  OF  CONVENTIONAL  RECEIVER 


1.  General  Error  Expressions 

As  we  pointed  out  previously,  what  we  mean  by  a  conventional  receiver 
is  one  which  is  designed  under  the  assumption  that  the  interference  is  additive 
white  Gaussian  noise.  Clearly,  when  this  assumption  is  correct,  the  resulting 
receiver  is  optimum.  Whenever  re ve liberation  is  present,  this  "conventional  " 
receiver  will  not  provide  the  optimum  processing.  Because  the  conventional  re¬ 
ceiver  frequently  is  far  easier  to  implement  than  the  optimum  receiver,  we  want 
to  find  out  how  far  from  optimum  the  conventional  receiver  is. 

In  this  case. 


T,(t)  =  ^  7(t) 

O 


(IV-48) 


Under  hypothesis  Hq, 


X  +  jy  =  J  ~  j 


(TV-49) 


-T 


-T 


T'  TT 

I  1  X  +  jy  I  ^  I  =  J*  f  dv  (f  (u)  q(v)  E  [rjj(u)  ^(v)] 


aV-50) 


-T  -T 


or 


Var[x] +Var  [y]  -  ^  f  du  j  dv  S*(u)  S^(v)  (u,"'  =  2o 

N  J  J  o 


av-51) 


-T  -T 


"ITten,  in  a  manner  identical  to  that  used  to  derive  Equations  IV-33, 
IV-38,  and  IV-39,  we  obtain 


c 


(IV-52) 
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Under  hypothesis  Hj, 


X  +  jy 


=  /  q*(t)  E  [rjw] 


dt 


-T 


/ 


-T 


tf  S^(t)ej-%t 

o 


2E 


N 


r  J- 


(IV-53) 


or 


?E 


N 


cos  = 


(IV- 54) 


and 


y  = 


2E 


Ik’ 

»N 


Sin 


o 


(IV- 55) 


Then 


H. 


(x.yl  =)  =  — 


exp  - 


(,.^cos=)  .(y-^sins)  | 


2rrC 


20. 


(IV- 56) 


Pt,  (x,y|s)  =-^  exp - 
T  2nc^ 


/2E  2E 


cos  3  +  ysin5 ) 


20, 


(IV-57) 
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,,» .  // 


Pu  (x.  y|  r)  dx  dy 


(lV-58) 


Changing  to  polar  coordinates  and  integrating,  we  have: 


Pd(=)  = 


CO 

i  /»/ 


/  /2E  \  2E 

Iz +^-^-2— ^Zcosce-?) 

dc  •  Zexp-I - - - - — - - - l(IV-59) 


The  answer  is  not  a  function  of  3,  so  we  have: 


i  / 


7?  + 


Pd  =  -  ZdZexp- 


(V) 


2E  Z 
r 


(IV-60) 


Just  as  in  the  case  of  the  optimum  receiver,  we  must  express  in 
terms  of  the  Q  function. 


/  2  3, 

(x  +  a  ) 


Q(a,3)  = 


I  (ax)  dx 
o 


(IV-61) 


Letting:  ^ 


=  X  in  Equation  IV-60,  we  obtain: 


X  dx  exp  - 


LVcJ 


O  NO, 


(IV-62) 
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ju 


- - -  «.U^ 

dt  dd  lit  uiw  wpviiituiAt 


&C  lO  «40V>1C»1 


t'r\ 

VW  WV'AAAIW 


SyV  2*2tiG- 


2 

d 

c 


tmean  of  output  1  H , 
_ I  1 

variance  of  output 


2 


(IV- 63a) 


Then,  we  may  write 


(IV- 63b) 


Using  Equations  IV-52  and  IV-63b,  we  may  write: 


P  =  Q  (d^.  -  2In  P  J 

d  c  F 


(IV- 63c) 


Thus,  the  quantity  d^  completely  characterizes  the  performance  of  the  conven¬ 
tional  receiver./  Fo:-  a  given  Pp,  we  may  use  Figure  9. 

We  now  v/ant  to  obtain  some  simpler  expressions  for  d^  for  the  sta¬ 
tionary  and  the  non- stationary  cases. 


2.  Stationary'  Case 

For  the  stationary,  infinite  interval  uase.  a  simple  expression  can  be 
obtained  in  terms  of  the  various  spectra. 

F'or  the  stationary  case.  Equation  IV-51  becomes 


0 


2 

C 


dv  S*(u)  S  (v)  R  (u-v) 
d  d  I 


/  This  was  pointed  out  by  C.  Boardman,  M.I.T. 


(IV- 64) 
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From  Equation  III-36, 


R  (u  -  v) 
n 

r 


y-£j  L?jJ  M  (u  -  V)  R  (u  -  V)  =  i  1  M  (u  -  V)  R  < 
2  Jj  I  2  av  r  I 


Defining, 


/ 


-jjut 


Sd(J^)  =  )  Sd(t)e  ^  'dt 


and 


=/ 


R(t)  e'-^^’^dt 


Then, 


du  S*  (u)  I  dv  e  (v)  (uj)  ^ 


Using  Equation  IV- 66,  we  obtain: 


f$ Iv I 


SJv)  SJi>) 
d  *  r 


Now,  from  Equation  III-31 


S  (uu)  =  N  +  S  (ID), 
r  On 

r 


where  S  (uj)  =  F  1  R(t  -t.)  1 

n  L  a  3  J 

r 
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-v)  (IV- 65) 


(IV-66) 


(IV-67) 


(IV-68) 


(IV- 69) 


(IV-70) 
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Using  ParsevaTs  rheorein,  we  have 


(IV-71) 


From  Equation  IV-63,  we  observe  that 


^  2E  N 
r  o 


a 

r  /  I  <"> 


(IV-72) 


3.  "valuation  of  for  Non-Stationary  Case 

conv  • 


From  Equations  IV- 51  and  III-31,  we  observe  that: 


2 

= 


(t  ,  t , )  dt  dt„  =  0  +0 
a  3  a  3  X  rc 


(IV-73) 


For  our  particular  problem,  the  desired  signal  is  an  atienuaicd  replica 
of  the  transmitted  envelope  Tit)  which  has  been  shifted  in  frequency  and  delayed 


in  tim.e. 


S^(t) 


2E^  f  (r  > 


(IV- 74) 


S^(t) 


2E  f(t-T  i  e 
r  d 


y^Tthur  fl.lLUtlc  Jnc. 


I  Then,  is; 


-jjj  t  +K,t 


(IV-75) 


From  Equation  111-23; 


t^)  =  (2E.)  ^  ifii'  Ja  (.)  M  (t^  -  tg:x)  T(t^  -  x)T  (tg  -  X)  dx 


(lV-76) 


We  recall  that 


"  +ji(t^-tg) 

1  (t  -t  :x)  =  I  p  e  dJ. 

q  I  -u:  Q I 


(IV-77) 


We  have; 


4E  E  i  e)  Izl^l 

_Lll_ - j  fj [dt^c-t,  dx  di.  T*  (t^  -  7(t3  -  T^)  r (t^  -  X) 


'^(t  -x)a(x)p  , 

t5  ‘X  I  X 


exp 


:iV-78) 
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_ -  Jfst- 


Rp-arrancrincr  rhp>  rf»rm«:  u)f»  hawo- 

- Q - Q - 


4EE  IeIIzPI 


t  r  2 


JdxJd 


ia(x)p 

ql 


/dt^T(t^-x)T*(ta-T^>exp  I  jt^Cx^-  x)j 
T  (tg  -  x)  f(t,  -  r^)  exp  j+jt^ 


(IV-79) 


First,  observe  that  the  fourth  integral  is  just  the  complex  conjugate  of 
the  third  integral.  Second,  observe  that  the  third  integral  is  identical  (except  for 
a  phase  shift)  to  the  two-dimensional  correlation  function  defined  in  Equation  III- 
24.  With  these  observations,  we  have: 


J dt„  T*  (t^  -  x)  f  (t.  -  exp  I  +  jt^  -  ii)  I 


8  (T  -  x;  X  -  x)  I  I  (T  i) 

d  d  d  d 


(IV- 80) 


Observe  that  we  normalized  T (t),  so  that  V  (0,  0)  =  1 . 

Equation  IV-80  is  just  a  definition  of  the  signal  ambiguity  function. 
See,  e.g..  Woodward  (Reference  1)  or  Siebert  (Reference  9). 

Next,  we  observe  that  the  term  a  (x)  p  ,  (j  :x)  is  just  the  joint  prob- 

q| 

ability  density  of  the  scatterers  in  delay  and  Doppler. 
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^  1  1  ^  ^4*^  «»«  ««  «^ 

VW  will  COiA  ClAiO  Cil\^  OV>a(.C\..  X  lil^  XUlA^i.X>.AiA 


s(x;x)  =  a(x)p  .  ^('^:x) 

X  I  X 

q  I 


(I.V-81) 


Substituting  Equations  IV- 80  and  IV- 81  into  Equation  IV-79,  we  have: 


4E  E  ^  E 1 1 Z  I 
t  r  2 


\ 


rc  '  d’  d' 


!N^ 


J*  dx  J*diJ  s  (x:ju)  V  (7^  -  x;  iw^  -  ui) 


(IV- 82) 


2 

Therefore,  for  die  conventional  receiver  aj.j;(Tjj;iDjj)can  be  expressed 
as  a  two-dimensional  convolution  of  the  reverbe’"'tion  scattering  function  and  the 
signal  ambiguity  function./ 

Combining  Equations  IV- 82  and  IV- 73,  we  have: 


,  2E 

/  = 
c  N 


^Ei  I  ZplE. 


1  + 


•  '  ’  t  rc 

- J  J  dx  du  s  (x:j>)  f  (T^  -  x;  -  a) 


(IV- 83) 


From  Eauation  IV- 63,  we  have: 


2E 


N 


1  + 


|e1  |z|^!e, 


(IV-84) 


N 


JJ dx  di>  s  (x:W;  f  (t^  -  x;  -  ia<) 


In  this  section,  we  have  derived  the  structure  of  the  optimum  receiver. 
We  observed  that  the  output  S/N  ratio,  d^,  provided  a  reasoaible  characteriza¬ 
tion  of  the  receiver  performance.  Frequently,  the  structure  oJ  the  optimum 


/  This  particular  form  is  not  new.  Westerfield  and  Stewart  (References  14  and 
iS)  obtain  a  similar  rtlation.  Green  (Reference  16)  describes  its  application 
to  radar  astronomy. 
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receiver  is  complex  or  requires  knowledge  that  might  not  be  available.  In  this 
case,  one  commonly  uses  a  convenaonai  receiver,  we  observed  that  the  per¬ 
formance  depended  on  tVk'O  parameters  and  o  To  keep  the  subsequent  work 
from  getting  immersed  in  details,  we  decided  to  use  the  single  parameter  d-  as 
a  basis  of  comparison.  The  principal  results  that  we  will  use  in  the  subsequent 
work  are  Equation  IV-44  (and  its  various  modified  forms;  e.g..  Equations  IV-46 
and  IV- 47)  for  the  optimum  receiver  and  Equations  IV-72  and  IV- 84  for  the  con¬ 
ventional  receiver. 

In  the  next  section,  we  will  derive  some  general  properties  relating  to 
signal  design  and  processing  in  the  presence  of  non- white  Gaussian  noise.  After 
deriving  these  properties  we  will  return  to  the  reverberation  problem  of  interest. 
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\7  CTi^TvTAI  oo  ^rir?0 'T*Tt?C 

V  •  ^  1^  i^CfOlV^JlN  £*  IW^I  L>i\  1 


I  In  the  preceding  sections,  we  have  formulated  the  detection  problem  in 

Ithe  presence  of  reverberation.  The  pertinent  results  can  be  summarized  briefly. 


The  optimum  receiver  formulated  -k  test  statistic 


1 

s  =  y*^ 


(V-l) 


where  q(t)  is  the  solution  to  the  integral  equation. 


1 


{V-2) 


The  performance  of  the  system  depends  on  a  quantity  d^,  where 


1 

d^  =  y*  q  (t)S^(t)dt 


(V-3) 


When  the  interference  has  a  “white”  spectrum,  i.e.. 


R  (t  ,tj  =  N  u  (t  -Q 
n.j,'  a'  s'  o  o  '  a  8' 


(V.4) 


^<'S>  =  -h  W 

o 


(V-5) 


>  2E 

1^  -  — i 

o  N 


(V-6) 


The  performance  depends  only  cn  the  received  energy  and  not  on  the 
signal  shape. 
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For  any  other  spectrum,  d  depends  on  the  signal  shape.  In  this  section. 


we  derive  several  properties  regarding  signal  design. 


A.  PROPERTY  1 

The  performance  of  tl^  optimum  ’’eceiver  is  minimized  by  choosing  the 
complex  envelope  of  the  signal,  Sjj(t  ),  equal  to  the  eigenfunction  of  the  noise 
with  the  largest  eigenvalue. 

The  proof  of  this  statement  is  as  follows.  For  any  threshold  Aq,  the 
performance  of  the  optimum  receiver  is  monotone  in  d^,  where 


i 

2  r 

i  =  I  a 

^  J 


it)  S^(t)  dt 


(V-7) 


X 


(V-8) 


Expand  q  (t^)  using  the  eigenfunctions  of  R  (t  ,t^) 

"t  ^ 


=1: 

1  =  1 


(V-9) 


where  the  0.  (t)  satisfy  the  integral  equation 
'  T 

-  T 


(V-10) 


Substituting  Equation  v-9  into  Equation  V-8  and  using  Equation  V-10, 


we  have: 


W  ’k“k^k<'> 


(V-Il) 


/  See  Reference  7 , 
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Similarly,  we  could  expand  using  the  eigenfunctions 


(V-12) 


Equations  V-11  and  V-12  imply 


'^k  ~  2 


(V-13) 


Now, 


/J  (  00  s  * 

He  ^0- 

“  I  1  =  1  0^  ■ 


<t)  E  s  0  (t)  dt 

j  =  l  1  1 


(V-14) 


Integrating  and  using  the  orthonormality  of  the  eigenfunctions,  we  have; 


(V-15) 


Now,  the  sum  of  the !  Sji  is  twice  the  energy  in  the  signal. 


setting; 


2E  =  E  I  Sjl 

i  =  l 


(V-16) 


2  2 

Denote  the  largest  eigenvalue  by  Oj^.  Then  clearly  d^  is  minimi  jied  by 


SlI  =  2E^ 


(V-17) 


Si  =  0 


i  t  L 


(V-18) 


B.  PROPERTY  2 

~  2 
If  S^(t)  is  chosen  to  give  the  minimum,  value  of  d^,  the  optimum  re¬ 
ceiver  is  identical  to  the  conventional  receiver. 


V'a>  =  '‘l  V 


(V-19) 


Arthur 
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Then,  Equation  V-8  tecomes; 

T 


-T 


F rcm  Equation  V  -9 ,  we  see  that : 


J2E 

q  (tg)  =  - 2^  0  (tg)  =  —  (tp) 

^  L 

which  is,  of  course,  the  conventional  receiver. 


(V-20) 


t 

j 


J 


(V-21) 


C.  PROPERTY  3 

^<w 

If  the  conventional  receiver  is  used  for  all  signals,  S^ft),  the  minimum 
performance  is  obtained  when  S^(t)  is  equal  to  the  eigenfunction  of  the  noise  with 
the  largest  eigenvalue. 

The  proof  is  analogous  to  that  of  Property  1 . 


D.  PROPERTY  4 

If  R„  (t  ,  t  )  is  a  positive-definite  correlation  function  corresponding 
‘  a  3 

to  a  process  with  finite  variance  and  the  0i.(t)  are  a  complete  orthonormal  set,  then 
there  is  no  unique  S^(t)  that  maximizes  dg  . 

The  positive-definiteness  implies  there  are  no  zero  eigenvalues  in  the 
expansion  of  the  noise .  Since  the  0j^(t)  are  a  CON  ^  set  there  is  an  infinite  number 
of  eigenvalues.  Thus,  there  is  no^mallest  one.  This  means  that  one  can  suc¬ 
cessively  increase  d^  by  choosing  Sj{t)  equal  to  eigenfunctions  with  successively 
smaller  eigenvalues . 


j  /  Complete  orthonormal . 

-  \ 

1 
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The  intuitive  meaning  of  this  statement  should  be  clear.  Suppose 
Rn^  (ti  ,ti)  =  0.  Then,  if  we  consider  the  typical  noise  spectrum  shown  in 
Figure  10,  we  see 


FIGURE  10  TYPICAL  NOISE  SPECTRUM 

that  the  eigenfunctions  are  of  the  form:^ 

(t)  =  k  cosab  t  -T<t<T  (V-22) 

n  n  n 

and  the  eigenvalues  are  of  the  form:: 

2  1 

0  =  - :j- 

"  4a  (1  +  b") 

n 

Where 

u2  u2  J-  J- 

12  3  n 

Thus, 

2  2  2 

a,  >a^  >-.->o  (V-25) 

12  n 

Now,  as  we  take  successively  smaller  eigenvalues,  the  corresponding 
eigenfunctions  are  cosines  of  successively  increasing  frequency.  Thus,  we  can 
make  the  system  perform  arbitrarily  well  by  transmitting  a  signal  of  arbitrarily 
high  frequency. 


(V-23) 


(V-24) 


/  See  pp.  99-101  of  Reference  7.  Note  that  there  are  also  eigenfunctions  of  the 
form  kj^  sinabj^t.  An  identical  argument  holds  for  these  eigenfunctions. 
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This  fact  was  obvious  from  the  shape  of  the  noise  spectrum.  The  lim¬ 
itation  here  is  a  practical  one  rather  than  a  mathematical  one . 

Observe  that  if  we  allow  processes  with  infinite  variance,  Property  4 
is  not  true . 

As  an  example,  consider  the  spectrum  shown  in  Figure  11 


FIGURE  11  AN  EXAMPLE  OF  NOISE  SPECTRUM  WHEN  PROCESSES 

HAVE  INFINITE  VARIANCE 


This  spectrum  has  a  smallest  eigenvalue, 

2  *r  2  ‘V  1 


N  o  N 

_ _ r_  2  _ _ £ 

min  2  *^1  2 


4a(l  +  b^) 


(V-26) 


However,  since  this  type  of  spectrum  will  not  arise  in  our  work. 
Property  4  will  be  valid . 

In  our  problem,  the  noise  correlation  function  is  determined  by  the 
transmitted  signal. 


R  (t.,t  )  =  R  (t  ,  t  )  +  R  (t  ,  t.,) 

Dt  1  2  »A  ‘  -  "r  ‘  - 


(V-27) 


R„  (t,.  t  >  =  2N 
"t  1  ^ 


J 

'o“o<‘ry 


a(x)M^  I  (t^-tg:x)S(t^-x)S  <tg-x)  dx 


-  T 


(V-28) 
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Thus,  as  one  changes  the  signal  shape,  one  changes  Rj^  (t  ,  t  ).  There- 

*1  *  ^ 


fore  the  problem  of  finding  the  worst  signal  shape  is  difficult  for  the  general  case. 
Moreover,  when  the  correlation  function  is  specified  by  Equation  V-28,  one  can  not 
always  satisfy  Equation  V-19.  We  will  confine  ourselves  temporarily  to  some 
simple  scattering  functions  and  signal  shapes.  Clearly,  the  reason  we  are  con¬ 
cerned  with  the  worst  signal  shape  is  that  it  tells  us  what  characteristics  we  want 
to  avoid  in  our  signal.  After  finding  the  minimum,  we  can  show  how  the  performance 
improves  as  we  move  away  from  the  minimum. 


In  Section  VI,  we  will  consider  a  stationary  reverberation  return.  In  Sec¬ 
tion  VII,  we  will  consider  a  non -stationary  return.  In  each  case,  we  will  evaluate 
the  performance  of  the  optimum  and  the  conventional  receiver  as  a  function  of  the 
signal  shape  and  scattering  function . 


,2hctfaiT  linlr  Jnf 


I  / 

i 

i 

i 

1 

VI.  GAUSSIAN  SIGNALS.  UNIFORM  REVERBERATION 
A.  ASSUMPTIONS  OF  THE  MODEL 


In  this  section,  we  apply  the  results  of  the  preceding  sections  to  a 
specific  situation.  From  our  discussion  in  Sections  III  and  IV,  we  observe  that 
to  evaluate  the  performance  of  the  conventional  receiver  in  a  uniform  reverbera¬ 
tion  environment  we  need  the  signal  shape  and  the  distribution  in  Doppler  of  the 
scatterers . 


i .  Signal  and  Receiver  Properties 

We  will  assume  that  the  velocit;/  of  each  scatter  is  a  zero-mean 
Gaussian  random  variable. 


Thus, 


„  rv\  -  ^  •'  B 

P„.  (X) - e 

D 


(VI- 1) 


where  B  is  the  rms  Doppler  shift. ^ 


To  use  Equation  ni-36,  we  require  the  characteristic  function: 


M  (jT)  =  e 
D 


(VI-2) 


We  will  assume  that  the  transmitted  signal  has  a  Gaussian  envelope  and 
linear  frequency  modulation. 


Thus, 


Re  [k  e 


S^(t)  “ 


The  complex  envelope  is: 

T(t)  =  ke 


2.2 

-at  -iDt 


(VI-3) 


(VI-4) 


/  Observe  that  B  depends  on  the  rms  scatterer  velocity  and  the  carrier  frequency 
of  the  transmitted  signal. 


I 

1 

I 

I 


1 

I 

! 

1 
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The  amplitude  and  instantaneous  frequency  are  shown  in  Figure  12.  We 
observe  that  the  pulse  duration  is  infinite.  This  is  an  idealization  which  makes 
the  analysis  appreciably  simpler. 


Since 


2E^  =  y  1  f(t)|  ^dt  j  ^  "2^" 


(VI-5) 


I 

—  I 


We  have: 


(tV)*  = 


(VI-6) 


To  evaluate  the  -performance,  we  need  the  Fourier  transform  of  the 
complex  envelooe. 


(VI-7) 


Completing  the  square: 


F(jJJ)  =  exp  - 


4(a+ jb)‘ 


lutegrating,  we  have: 


(VI-8) 


F(jr)  =  2E; 


(2(a+jb))' 


4(a+jb) 


=  2E^  (2^)* 


2(a^+b^ 


exp  - 


2  2 
4(a  +b  ) 


2 

uu  -  jb  1 .  -lb 

- I 

4(a  +b  )  ^ 

(VI-9) 
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ENVELOPE 


INSTANTANEOUS  FREQUENCY 


FIGURE  12  MAGNITUDE  AND  PHASE  OF  COMPLEX  EN’VELOPE 


If  we  define: 


2  2 
(a  +  b  ) 

a 


(VI- 10) 


we  may  write: 


F(jj))  =  2E^  (2n)* 


1 


i  I  ■ 

Ci-fy 


Ui 

"4 


r-.p 


b  1  -1  b 

2-2--^“"  I 


4(a  +b  ) 


* 


(VI-11) 


The  magnitude  squared  is: 


F(ju;)|^  =  4E^  -  (2n)^ 


2  A 


1  'i/ 

-jexp- 


2A 


-OP 

-  2^t  ‘  rrj 


(VI- 12) 


As  a  simple  check  on  the  constant,  we  observe  ihac 

2  dui 


/ 


|F(j>-)| 


,  =  2E 

2n  t 


(IV- 13) 


Now,  the  correlation  function  R^('i‘)  is  simply  the  inverse  transform: 


9  2 

(^)  =  J  2E^  exp  -  Iy  exp  -  ym 


duu  _  AT 

2n  t  ^  2 


(VI- 14) 


From  Equation  IV-65,  we  have  (for  the  stationary  case): 


I 

av 


R  D 


(VI- 15) 


Substituting  Equations  VI-1  and  VI-14  into  Equation  VI-15,  we  have: 

R  (T)  =  E  •  I  exp -4(B^+A)t^  (VI-16) 

n  t  av  ^  2  '  '  ' 

R 

Let 


2  „2 
Y  =  B  +  A 


(VI- 17) 
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t 


Sf- 

•t 

4’. 


Than, 


2  2 


R  ('^  =  E  •  I  exp  - 
n  t  av 


(VI- 18) 


and  the  reverberation  spectrum  is; 


2E  ♦  I 


-  /.X  t  av  v2n  u; 

<“>  =  —2 - 7-  “P  -  —2 

R  2  Y 


(VI- 19) 


Using  the  conventional  matched  filter,  we  have; 


12  (“"'V 

=  2E  — X—  exp  - 


(VI -201 


2.  Evaluation  of  Conventional  Receiver 

To  evaluate  the  conventional  receiver,  we  sirnply  substitute  Equations 
\T-19  and  VI-20  into  Equation  IV-72.  ITiis  gives; 

2 


d^  = 
c 


2E 

“n 


r  1  f  ../2^ 

o  K  J 


(a  -  a:  )  2E  •  I 

D  t  av  72 n 


exp 


2v 

(\T-21) 


JL.  ^ 

2  2tt 


or 


2E 


d“  = 


N 


E 

!-{-  -J 


+  3D 

T^av  1  f  ( 


2A 


+ 


'JU 


dir 


2v 


-1 


(VI-22) 
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Collecting  terms  and  completing  the  square,  we  have; 


2E 


N 


d  = 


^/^T^ 

“i 


fJl^ip-poA 


1  +  2E  -  I  exp  - 

t  av  Y  2Y^ 


da; 


where 


2  2 
Y  =  B  +  A 


Evaluating  the  denominator,  we  have: 


den.  =  1 + 


E  •  I 
t  av 

N  A^Y 


w 

/ 


exp  - 


Tr~ 


exp  - 


U) 


2y' 


Letting 


c  - 


N 


0/2 


and  collecting  terms  and  completing  the  square,  we  have: 


den .  =  1  + 


E,s  ; 

2N  A^y  J 


exp  - 


B^4-2A 
2A(B^+  A) 


B 


b^+a1 
B^+2  AJ 


(VI -23) 


(VI-24) 


duj  (VI-25) 
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duu 


•  exp  - 


2(B  +2A) 


(VI-27) 


Integrating,  we  have: 

Tru  A^  (B^+ A)^ 

den.  =  I  -r  — T -  - 2 - Z -  ®^P  *  7 

2A^v  (B  +2A)^  2(B  +2A) 


D 


(VI-28) 
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This  expression  reduces  to: 


H  6  n 

den.  =1+  — ^ - T-  exp - z -  (Vl-29) 

2(B  +2Ar  2(B  +2A) 


Now  consider  the  behavior  as  a  function  of  uUj^,  B,  and  A  . 
For  B  0,  this  can  be  rewritten  as: 


den. 


1  + 


E  ByTfr 
t 

2B 


(VI-30) 


We  observe  that  the  ratios  of  importance  are: 


'\)  :  The  ratio  of  the  target  Doppler  shift  to  the  rms 

B  reverberation  Doppler  shift 


and 

VTa^ 

B 


:  The  ratio  of  the  effective  bandwidth  of  ti.e  signal 
to  the  rms  reverberation  Doppler  shift . 


To  maximize  d^,  we  want  to  minimize  the  second  term  in  Equation  VI-30. 
The  first  coefficient  is  a  function  of  the  environment  and  the  transmitted  energy. 
Considering  this  constant,  we  want  to  study  the  behavior  of  the  function: 


The  curve  is  shown  in 


Figure  13. 


/ 


/  One  observes  the  similarity  between  Figure  13  and  Figure  1  of  Reference  14. 
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FIGURE  13  PLOT  OF  f 
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TTie  minimum  value  of  f 


fr-4^) 


occurs  when 


»d<b 


(VI-32) 


=  0 


a,p<B 


At  the  minimum: 


{?■ 


J2l? 


'^T  /  2  /  \  ^ 

=  1.65—  =  1-65  + 


(VI -33) 


Clearly  the  importance  of  the  term  given  by  Equation  71-3 1  depends  on 
the  value  of  its  coefficient  in  Equation  VI -30. 

To  study  this  effect,  we  substitute  Equation  VI -30  into  Equation  VI -22: 


d^  = 
nc 


1+D 


2/  2o’ 

2B  (l+-j 

\  B  , 


(VI-34) 


where 


E  BjTn 
t 


(VI -35) 


Physically  this  represents  the  reverberation-to-ambient-noise  level  in 
the  reverberation  bandwiddi . 


E  I  .  N 


Recall  that  8  = 


.  Thus 


w  A  I 


8 

V  2rr 


The  denominator  is  simply  the  noise  power  out  of  a  filter  with  a 
Gaussian  spectirum  when  the  input  is  white  noise  of  spectral  height  Nq/2-  From 
Equation  HI -32,  we  observe  fiiat  the  numerator  is  just  the  total  received  power. 
Thus,  D  is  die  reverbcrarion-to-ambient-noise  ratio. 
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I 

4r 


I 

I 
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The  physical  meaning  of  should  be  clear.  In  the  absence  of  rever¬ 
beration  it  equals  1.0.  A  decrease  to  a  value  of  less  than  1 .0  represents  the  loss 
in  detectability  due  *  reverberation . 

We  have  plotted  following  para«T^eters; 


D  Physical  Meaning  Figu 

0.3  reverberation  < ambient  noise  14 

1.0  reverberation  =  ambient  noise  15 

10.0  rev. /ambient  noise  =  10  db  16 

100.0  rev. /ambient  noise  =  20  db  17 


The  parameters  on  the  curves  are  This  is  the  ratio  of  the  target 

velocity  to  the  reverberation  Doppler . 

1 

The  horizontal  axis  is  L^/B.  This  is  the  ratio  of  die  signal  bandwidth 
to  the  reverberation  Doppler. 

Several  observations  may  be  made  with  respect  to  this  class  of  signals: 

a.  For  zero  target  velocities,  we  have  monotone  improvement  as  the 
bandwidth  increases. 

b.  For  non-zero  target  velocities,  one  can  use  either  very  small  or 
very  large  bandwiddi  signals ,  The  point  of  the  exact  minimum  is  a  function  of 
Wd/b  given  by  Equation  VI-33. 

c.  For  small  L,  the  non-zero  target  velocities,  the  improvement  in¬ 
creases  rapidly.  However,  for  large  A,  all  targets  behave  the  same.^ 

Now  we  consider  the  optimum  receiver. 


/  This  statement  is,  of  course,  for  the  class  of  signals  that  we  are  considering 
(i.e. ,  linear  FM  and  Gaussian  envelopes).  It  should  be  emphasized  diat  this 
statement  is  not  true  for  all  signals. 
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FIGURE  14  LOSS  IN  DETECT AIIILITY  DUE  TO  REVERBERATION  (D-().3) 


FIGURE  16  LOSS  IN  DETECTABILITY  DUE  TO  REVERBERATION  (D=  10.0) 


t 


I 


B.  PERFORMANCE  OF  OPTIMUM  RECEIVER 


We  will  use  die  same  assumptions  as  in  the  preceding  section .  Now, 
from  Equation  IV-46-  we  have 


(VI -3  6) 


Substituting  Equations  VI- 19  and  VI -20  into  VI -36,  we  have: 


This  can  be  rewritten  as: 


(VI -37) 


(VI-38) 


For  arbitrary  parameter  values,  one  cannot  obtain  a  closed  form 
solution.  We  will  first  consider  the  case  where  the  reverberation  is  small. 


1 .  Perturbation  Solution  (Low  Reverberation  Levels) 


Consider  the  case  when 


E3 

t 


Jlv 


<<  1 


(VI-39) 


In  this,  the  denominator  can  be  expanded  in  a  series  which  is  absolutely  con¬ 
vergent  . 
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a 

3k 


I 

I 


Then, 


o  2E 
.2  _  t 

N 

o 


V2n 


/ 


exp 


2A 


f  I 

..  I 

I 

t 

i 


y2iT^  (  op^  L^t'®  -/itiN  l  2j)' 
(_._jexp  --^M—  — jexpj-- 


+  •••  (-1) 


2E 


[-2 - -)  “P  j- 


nyb 


ly 


The  first  term  is  just  tt—  .  the  conventional  filter  result. 

N 

o 


The  second  term  is 


o  2E 


E  -B 

t 


du) 


(VI-40) 


V2tt 


u 

(id  -  U)_  )  2 

'  D  ou 

2/i  .,2 

J 

go 

2  Y 

du) 

(VI-41) 


Completing  the  square  and  integrating,  we  have: 
-2E 


1  •''2n  ^4.,  ^ 

-Z—r  •  -2 - ^  •  —j - i-‘/2TTexp - 2 - 

o  V2ttA-  ^  ^  (B+2tr  2(B  +  2A) 


I 

I 

I 

I 


■y 


'2rr 


'JU, 


D 


“2 - i"  ■  - 2 - 

(B  +2A)®  2(B  +2  A) 


Similarly,  the  (n+l)st  term  is: 


exp  - 


(tU-UUj^' 

TZ 


nu) 

2N'^ 


(VI-42) 


duu  (VI-43) 
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Completing  the  square  and  integrating,  we  have: 


.  2E  /E  -8  nr~\^ 

o 


Then, 


(-1)"  Y 


(B  -(iH-1)^)' 


D  n 


exp 


B  +(n+l)  L 
(VI -44) 


2E  ®  /E  8  rr-\ 


exp  - 


(B  +(iH-l)A)‘ 


2  2 

B  +(tH-l)  L 


(VI-45) 


Using  the  same  type  normalization  as  previously  we  can  write:  (for  Bj^O) 
2E_  (  /E  S  yTS'X  .  .  /m  . 


2E  1  1  ifoV 


(E  B  y2n\ 

(..A) 


exp  - 


(t/ttsii 


(VI -46) 


We  observe  that  the  first  two  terms  of  the  series  are  identical  to  a  series 
expansion  of  d^ 

2E  E  3  yyr  2 

“c"  TT^  1-  _  exp  (Vl-47) 

O  2  Y 


which  is  the  signal-to-noise  ratio  at  the  output  of  the  conventional  filter. 
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Several  observations  may  be  made  by  examining  Equation  VI -46: 


a.  As 


-  ®,  the  performance  of  the  conventional  and  that  of  the 


optimum  filter  will  approach  each  other.  Physically,  this  means  that  a  signal 
whose  bandwidth  is  much  greater  than  the  reverberation  Doppler  tends  to  elimin¬ 
ate  the  reverberation  disturbance. 


b.  As 


'  - ,  die  performance  of  the  conventional  and  that  of  the 


optimum  filter  will  approach  each  other.  Physically,  this  means  that  if  the  target 
Doppler  is  much  greater  than  the  re  verberation  Doppler,  the  effect  of  the  rever¬ 
beration  is  reduced. 


c.  As 


E  B 

t  v2tt 


0  ,  the  performance  of  the  conventional  and  that  of 


the  optimum  filter  will  approach  each  other.  This  is  just  die  obvious  fact  that  as 
die  reverberation  strength  goes  to  zero,  the  optimum  filter  becomes  the  conven¬ 
tional  filter. 


These  observations  really  indicate  when  the  effects  of  reverberation  are 
not  important .  Now  consider  the  case  for  arbitrary  reverberation  levels . 


I 

I! 

I I 


2.  Arbitrary  Reverberation  Levels 


Now  we  consider  the  general  case.  The  normalized  d^pj.  is: 


2E  I  E  B  /T — 

« «  A  '  A  %  • 


exp  - 


(VI-48) 


UU 

Letting  x  =  — and  rewriting,  we  have: 


no 


exp  J -  2  I  X  - 


_L_/ - L_ 

J  E^e  ^ 


A^/Bj 


1  X 

2 


(VI -49) 
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We  ob.serve  that  the  quantities  of  importance  are; 

ratio  of  target  Doppler  to  reverberation  Doppler 

ratio  of  effective  signal  bandwidth  to  reverberation 
Doppler 

A 

=  D  ;  ratio  of  the  total  reverberation  power  level  to 
the  noise  power  level  in  the  reverberation 
bandwidth.-*^ 

The  integral  in  Equation  VI-49  was  evaluated  numerically  for  several 
cases.  The  results  are  shown  in  Figures  18  and  19.  We  observe  that  the  curves 
have  the  same  general  characteristics  as  the  curves  for  the  conventional  receiver . 
For  D  =  0.3  and  1.0,  the  quantitative  difference  in  the  performance  is  insignifi¬ 
cant  and  the  curves  are  the  same  as  those  in  Figures  14  and  15.  However,  for 
D  =  10  and  100  (i.e. ,  reverberation/ambient  noise  ratios  of  +10  and  +20  db), there 
are  appreciable  differences  in  some  cases . 

Let  us  examine  the  effects  of  these  differences .  To  demonstrate  the 
method  of  conqjarison,  consider  a  specific  example,  in  which  the  parameters  are; 

a.  =  5,0 

b.  A^'/B  =  1.0 

c.  D  =  100.0 

For  this  set  of  parameter  values; 

d^  =  0.762 
no 

and 

d^  =  0.528 
nc 


(VI -50) 

(VI-51) 


a .  /B  ; 

b.  a^/b  : 


c. 


E  S  J2t\ 

t 

2B 


/  Recall  that  9  =  - 

0/2 


^t'  /  ^  3 

Thus,  D= 

2  •  2  /2n 


The  denominator  is  simply  the  noise  power  out  of  a  filter  with  a  Gaussian  spec¬ 
trum  Mv4ien  the  input  is  white  noise  of  spectral  height  ^Qf2'  From  Equation  III -3  2, 
we  observe  that  the  numerator  is  just  the  total  received  power.  Thus,  D  is  the 
reverberation-to-ambient-noise  ratio. 
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FIGURE  19  OPTIMUM  RECEIVER  PERFORMANC 


One’s  first  reaction  is  that  since  the  difference  is  about  1 .5  db,  tlie  value  of  an 

rkr%»imtirrt  fllfov*  ic  ni  to  “sH /\na  Klo  r*  rHic  r^om  r«a  i*i  crtn  ic  crvrrtot&fHaf  rt'itc- 

leading.  Let  us  assume  that  we  recniire  a  given  d^  to  obtain  a  desired  performance 
level.  Denote  this  required  d^  by  d^  .  Assume  that,  with  the  above  param¬ 
eters,  the  optimum  system  provides  the  required  d^  .  Thus, 


2E  - 

-rr^^  = 

N  r 

o 


(VI-52) 


Now  we  want  to  find  how  much  we  must  increase  the  transmitted  signal 
energy  to  achieve  the  same  d^  with  conventional  processing.  If  we  increase  the 
transmitted  energy  by  a  factor  k,  we  have; 


E  =  k  E 
rc  ro 


(VI-53) 


E  =  k  E 
tc  to 


(VI -54) 


From  Equations  VI -34  and  VT-48,  w'C  have: 


no  N 


1  4-  Dk 


For  the  above  case; 


L  ^ 


(-7 


I  (V* 

^^^-0  7  -  ^ 
*'  fl  +  2^ 


(VI-55) 


1  -  d  •  D 
no 


\  ^ 


exp  - 


1  IVBf 

2  /  2> 
(1  +  2  VB 


(VI -56) 


(If  it  is  impossible  to  achieve  the  required  d^  by  increasing  the  transmitted  energy. 
Equation  VI-54  may  have  a  negative  solution.  Clearly,  k  must  be  positive  to  be 
meaningful.) 


>• 
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tivaiuacing,  we  obtain: 

k  =  2.45  (+3.9  db)  (VI-57) 

This  increase  in  energy  required  when  a  conventional  receiver  is  used 
instead  of  an  optimum  receiver  gives  a  more  accurate  picture  of  the  cost  of  non- 
optimum  processing.  The  difference  in  the  two  results  occurs  because  as  the 
transmitted  energy  increases,  the  reverberation  return  also  increases.  Thus, 
for  hig^  reverberation  levels,  one  cannot  combat  the  reverberation  by  raising 
die  energy  level . 

The  conclusions  to  be  drawn  from  tiiis  section  are  two-fold: 

a.  The  most  important  step  in  combating  reverberation  is  proper  signal 
design.  This  signal  design  can  be  accomplished  using  a  conventional  receiver. 

As  pointed  out  in  Section  V,  one  will  not  achieve  a  unique  maximum.  The  allowable 
range  of  signal  parameters  will  be  governed  by  such  factors  as  available  bandwidth, 
pulse  duration,  sound  path  stability,  and  other  system  constraints.  The  choice  of 
signal  will  also  be  governed  by  what  range  of  target  velocities  is  of  interest. 

b.  The  second  step  is  to  consider  how  much  one  would  gain  by  using  an 
optimum  receiver  instead  of  a  conventional  receiver.  One  must  examine  the 
specific  situation  to  see  if  the  added  complexity  is  warranted . 

Now  we  want  to  conduct  a  similar  analysis  for  the  non-stationary  case. 
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MnM_f  IMICODKA  ODVCDDDO  AT^rkKI 

WA«4A  V..#*\A«A  »  *.  4  « 


In  this  section,  we  consider  the  performs  ce  of  conventional  and  opti¬ 
mum  receivers  in  the  presence  of  non-uniform  reverbe ration »  As  discussed  in 
Section  III,  if  the  distribution  in  range  of  the  scatterers  is  non-uniform,  then  the 
reverberation  return  will  be  a  sample  function  from  a  non- stations iry  Gaussian 
process.  We  consider  first  the  conventional  filter  performance.  In  this  section, 
we  restrict  ourselves  to  pulses  with  Gaussian  envelopes  and  linear  FM. 


A.  CONVENTIONAL  FILTER 


As  pointed  out  in  Section  IV,  to  evaluate  the  conventional  receiver,  we 
require  the  ambiguity  function  of  the  signal  and  the  scattering  of  the  reverbera¬ 
tion.  First,  evaluate  the  ambiguity  function. 

From  Equation  III-24,  we  have: 


4® 

(T,<r)  =  J  7  (t -j)  7*  (  +j) 


-  jxt  j 
e  dt 


(VII- 1) 


For  a  Gaussian  pulse. 


2E,  f(t)  = 


/saE^V  i  2  1 

1  I  exp  I  -  (a  +  jb)  t  I 


First,  let  b  =  0. 


Then,  substituting  Equation  VII-2  into  VII- 1,  we  have: 


9(t,x)  = 


(-J  / 


exp  I  -a  (t-^)  -a(t+^)  -  jJJt  |  dt 


(VII-2) 


(VII- 3) 
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Completing  the  square  and  integratinj: 


9(T,a>)  =  2E^  exp 


3  2 

aT  uu 

'2~‘  8a 


(VII- 4) 


Then, 


¥(T,a))  =  |e(T,uj)|  ^  =  4E^  exp  j-  ^  | 


(VU-5) 


Normalizing, 


„  .  .  A  Y(T,uj)  )  2  lu 


(VU-6) 


A  sketch  of  the  ambiguity  function  is  shown  in  Figure  20. 
For  b  5^  0,  we  obtain:/ 


9(T,a,)  =  4E^  exp  j  -  2^  - 


(vn-7) 


X  ,2  '  2bT) 

'i  fT ,  uu)  =  exp  ( -  aT  -  - - - - 

n  I  4a 


(VII- 8) 


The  equal  amplitude  lines  for  the  linear  FM  cases  are  shown  in  Fig¬ 


ure  21. 


/  One  can  either  evaluate  the  integral  or  use  Theorem  4,  Siebert  (Reference  17). 


3tiI|ut  2l.%UtIf,lnc. 


FIGURii  21  AMBIGUITY  FUNCTION,  LLNEAR  FM  PULSE 
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We  assume  that  the  scattering  density  is  a  skewed  Gaussian  density. 


s  (x.  Z)  = 


2n\  1  -  p  LB 


2  2  2  2 
1  B  X  -  2BLp  xZ  +  L  Z 

exp  -  X  - - 

^  2  2  2 
B  L  (1  -  P  ) 


(VII-9) 


The  simplest  case  is  p  =  0  (no  skew). 


Then, 


(vn-io) 


Substituting  Equations  VII-6  and  VII- 10  into  Equation  lV-84,  we  have: 


c  , 


jb{|z|^(e 

- - J  j  dxdii  s(x,uu)  ^  (T^ -  x;  -u;) 


-  2E  IeIIzI^Ie 

_  r  ,  2 

^^c  N  P  '  ■  N 

o  I  o 


!(uu 

-a(Td-T)  - 


4a  ) 


(Vn-11) 
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^  .  ..U.-.,..,. .  .  ■■  ■>■,,•■.  .  . . .  ■  -  - — - - - - - - - -  n 


ii 


£ 

i 


Evaluating,  we  have: 


d  ('f  ./A  .) 

c  d  d 


2E 
_ 1 

N 


1  + 


jbI  IzI^Ie^ 


N 


exp 


]_ 

2 


r  2 

‘d 


LU 


,2  1 


r,2  ^ 

B  +  2a 


-1 


(VII-.12) 


First  look  at  the  behavior  of  (0, 0)  as  a  function  of  a , 


The  worst  case  is  when 


_  B_  A 
^  "  2L  ®x- 


(VII-13) 


Then, 


d"  (0. 0) 
c 


2E  I  -e{|z|^}e 

_J[  1  +  2 _ '_t  _ 

N  )  N  BL+1 

o  I  o 


-1 


1 


(VII -14) 


A  convenient  way  to  examine  the  behavior  of  d^^^^  (0, 0)  for  other  a 
is  to  define: 


“  =''2BL= 


and  examine  In 


d"  (0, 0) 

conv 

2E  /N 
r  o 


=  In  d^  (0,0) 
nc 


(VII- 15) 


(VII -16) 
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Thus, 


(0,  0)  =  + 


iellzl^lEj 

N 


1  RI  ■■ 

(1  +  kBL)''  (l+-j^) 


(VII- 17) 


-  BL  ^ 

In  Figure  22,  we  plot  (1  +  JcBL)^  ^  "k*^  ^  function  of  k  for 

BL  =  1.  Since  the  function  is  symmetric  (on  a  In  scale)  about  k  =  1,  we  need  to 
plot  only  k  ^  1 .  CXiserve  that  the  nearer  to  one  the  value  of  this  function,  the 
better  the  system  performance. 

In  Figures  23,  24,  25,  and  26,  /  we  plot  d^^  (0, 0)  for  various  values  of 

^e{|zI^}e 

- - r; - -  ,  BL,  and  k. 

N 

o 


For  non- zero  and 


uu 


D’ 


we  want  to  choose  a  to  make  the  term 


as  small  as  possible. 

As  before,  there  is  no  unique  minimum.  There  is  a  worst  value  of  a. 
This  worst  value  is  a  function  of  3,  L,uu£j,  and  Tp.  There  does  not  seem  to  be  a 
simple  analytic  expression  for  the  worst  value. 

Similarly  one  can  use  the  ejqiression  in  Equation  VII- 8  for  the  am¬ 
biguity  function  of  a  Gaussian  pulse  vddi  linear  FM. 


/  Tlie  optimum  filter  result  is  also  plotted.  This  will  be  derived  in  Section 
VII- B  and  discussed  at  that  time. 
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28 


FIGURi-  '3  PERFORMANCE  OF  CONVENTIONAL 
OPTIMUM  RECEIVERS  — BL  =  0.10 


B.  OPTIMUM  FILTER 


rived  in  Sea,  “  *'-tegrel  equation  de- 


QD 

:V  "  /  '*r  V'3>^<'8>  *B 


(VU-18) 


where 


®r<'a-'3> 


^o‘<‘a-V"«n  <'a-'3> 

r 


(VIM9) 


’'n  <‘a-'s> 
r 


I  2^UZ\^] 


2i  *  /t  +tA 

’  f  +M-^) 

“j®<'^-'a-y*'<‘s-V“'>^  V  2  1^ 


(vn-20) 


iTEdZl^}  2E 
,,tj  =4^— - - - L 


R  ft .  1 1 = -  —  '  t  1  9  r  9 


(B^  +  a)  +  —2 


P-i) 


r,„2 


(B  +a)  - 


4fl.^-L 

\  4a 


<vn-2i) 


I 


Tc  solve  the  integral  ennatioTi;  we  wiii  evpanH  the  kernel  in  3  suitable 
bi-orthonoriml  expansion. 

A  suitable  series  can  be  obtained  from  Mehler's  expansion,  /  which  is: 


2  2, 

X  -i-y  -2pxj 

2(1-0^) 


,  r  2  2*]  »  H  (x)  H 

1  X  y  n  n  n 

=  ^  exp  -  —  -  ^  ^  - T  •  — 

^'L  ^  "  (n!)"  (n!) 


I  0  I  <1 


where  Hjj(x)  is  the  n^  Hermite  function. 


Since 


y[~^ 


OD 

r/ 


H  (x)  H  (x)  exp 
m  n 


dx  =  n! 


n  =  m 


otherwise 


A  suitable  set  of  orthonormal  functions  are: 


0  (x)  =  - X - r  H  (x>  exp 

(2*t)*  (n!)" 


From  Equation  VII-22,  we  obtain  the  desired  expansion: 


1  x^(l+p  )  +  y^(l  +  a  )  -  -^xy  n  ^  ^  , 

- _  exp - i - ^ - ^  =  2^  -  0„(x)  0  (y) 

4(1-p^  n=U  "  " 


We  observe  that  all  of  the  eigenvalues  of  the  Gaussian  kernel  are 
specified  by  a  single  number  P . 

/  See  Reference  18.  Our  approach  here  is  related  to  that  in  Reference  19. 


(\ai-22) 


(Vn-23) 


(Vn-24) 
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Now  we  want  to  expand  R  (t  ,t„)  as  given  by  Equation  VII-2i  into  a 
^  n  a.  c 
r 

series  of  the  form  of  Equation  VII -25. 

Introduce  two  new  variables; 


2  2 
r  =  R  + 


(VII-26) 


(VII -27) 


Then, 


R  (t  .t  ) 
n  a  3 
r 


1  1 

2 


^E 


4DE^  exp 


I  \tl(C^  +  D^)  +  tl(C^  +  D^) 


-  2t^tg(C^-D^ 


(VII-28) 


Now  let  X  =  ta,  y  =  t„  a  and  solve  for  o  and  P  such  that  Equation 
VII-28  is  identical  to  Equation  VII-22. 

We  obtain 


C  -  D 
C  +  D 


(VII- 29) 


(VII- 30) 


After  a  little  algebra,  we  obtain: 


R  (t  .tj  = 
n  a  3 
r 


iEdZl^ir  8aE^ 


n(l-^4aL  )J  (C+D)' k 


V  [(CJD)] 

^  [(C+D)J  k^E 


(\TI-31) 
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where 


fv.  u. 


(VII-32) 


VII-31: 


Including  the  white  noise  term,  we  obtain  from  Equations  VII- 19  and 


(VII-33) 


where 


iE(lzl^)r  8aE^' 


”(l  +  4aL^)  (C+D)^ 


(VII-34) 


Now  we  define  an  inverse  to  R  (t  , 

r  a  3 


OD 

/«r<'a-'9>'"r<'3’V  “'s  = 


^  (t  -  t  ) 
a  Y 


(VII-35) 


It  is  easy  to  verify  that: 


r  (t  ,tj 
r  a  3 


rT  '('a-'s> 

O 


4  2  — 


o  n=l  N  +K  0 
or 


(VII- 36) 
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uLMKnu  me  coeiucienr  in  tne  sum: 
n 


K  . 

C  ^ - £ - 

"  N  +  K 
o  r 


N 

n  o 

^  *r 

r 


This  gives: 


Multiply  both  sides  of  Equation  VII-18  by  T  (t  .t^)  and  integrate. 

r  Cl  r 


00 

/  Sd('a>"r<'a '  / 1  "r  W  V^<V  ‘‘S 


which  reduces  to: 


Substituting  Equation  VII-36  into  VII'*39,  we  obtain: 


00  ® 

q(t  )  =  4-  y*  C  ^(1)1:  (tj  S  ftj  dt. 

N  d  a'  N  /L/  n  n'a'  I  n  ='  d  z  c 

o  o  „_n  J 


00 

=  /  ’n<‘3>  Sd<'3>  “'s 


(VII-37) 


(VII-38) 


(VII-39) 


(VII- 40) 
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Then. 


q(t  )  =  T7-  S. 

*  a'  N  d 


S .  (t  )  -  ^  2  C  r  : 

M  u  a  N  n  n  n 

o  o  n-1 


•  (t  ) 


(VII -42) 


which  specifies  the  optimum  detector.  It  is  simply  a  band  pass  filter  whose 
complex  impulse  response  is  matched  to  q(t). 

To  evaluate  the  {jerformance.  we  recall  from  Equation  IV- 44  that; 


C  =  S.(t)q*(t) 

o  Ida  a 


(VII- 43) 


Therefore,  using  Equation  VII- 42,  we  obtain; 


-  ± 
^0  ■  N 


00  00 

/s,(t  )S*(t  )dt  -  Y'  C  r  I  f  (t  )  S^(t  ) 
d  a  d^a  a  nn  J  n  a  d  a 


(VII -44) 


which  reduces  to: 


'o  =  TT  -  if  2  cjrj 


o  n=l 


(VII- 45) 


Now,  the  are  a  CON  set,  so  we  may  write: 


Sd<‘a) 


(VII -46) 
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which  implies: 


S  (t  )  S*(t)dt 
d  d  a  a 


(VII- 47) 


Since  the  eigenvalues  are  monotones  decreasing  with  increasing  n,  it 
follows  from  Property  1  of  Section  V  that  d^p^  is  a  minimum  when 


r  =  0 

n 


(VII-48) 


S,(t  )  =  r  :  (t  ) 

d  1  o  o  3- 


(VII- 49) 


We  see  that  this  implies: 


exp  <  -  a^t  >  =  exp 


2  V 

16a  (B  -i-a  ) 
w  w  J 


4a  L 
w 


(VII -50) 


w  2L 


(VII-51) 


From  Properties  2  and  3  of  Section  V,  we  know  that  Equation  VII-51 
must  be  identical  to  Equation  VII- 13.  In  other  words,  the  worst  signal  is  identi¬ 
cal  for  the  conventional  and  optimum  receiver. 

To  evaluate  the  error  for  the  worst  case,  we  must  compute  Cq  and  r^. 

From  Equations  \1I- 37  and  VII-47,  it  is  clear  that: 


r  ^  J~2E 
o  f  r 


(VII-52) 


^Tt!)ur  m.Hittlt.inc. 


c 

0 


1 


(VII- 53) 


C 

o 


N 


1  + 


o 


K 


and  from  Equation  VII-34 


K 

r 


1. 

2 


E  {  I  Z  Pi  E 

t 

BL+1 


(VII-54) 


Then, 


ow 


2E 


N 


^E 


1  + 


N  BL+1 
o 


(VII- 55) 


(VII -56) 


Equation  VII-56  gives  the  same  result  as  Equation  MI- 14  since  the 
filters  are  identical. 

For  the  general  case,  the  optimum  filter  and  the  conventional  filter 
will  be  different.  To  evaluate  the  error,  we  must  sum  the  series  given  by  Equa¬ 
tion  MI- 45. 

Writing  a  =  ka^^,  we  have: 


K 

r 


i  E[1z|^]E^ 


i  1  1 

(k  +  2BL)’  (1  +  k2BL)?  +  k^ 


1 

? 


(MI-57) 
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and 


I 


I 

i 

I 

I 

I 

I 

I 

I 

I 

I 


i 


(k  +  2BL)^  (1  +  k2BL)^  - 

-  ■  '  ii 

(K  +  2BL)^  (1  +  k2BL)’  +  k^ 


and 


c  =  f  , 

n  J  n 


(t)  S^(t)  dt 


Now, 


s,(t) 


■(^) 


-at 


and 


H  (t) 
n 


where 


H_^[t[4CD]^]e‘'  “ 

{hY 


i  [2BL+k] 


_2  A  B  1 

D  =  :rr  r 


2L  (1  +  2BLr) 


95 


(VII -58) 


(ViI-59) 


(\ai-60) 

i 


(VlI-61) 


3 


I 


(VII -62) 


(VII- 63) 


t 


-V  "--<5-'  , 


Then, 


(16E  ;  a 


n 


&E  )  r  r  1 

rr/TTi  J".[‘'‘“>T'' 


(n!)^  (2tt)" 


( 


4CD  I 


-t  CD^-at  (vn-64) 


Letting  x  ~  t  [4CD]"',  this  reduces  to: 


r  = 
n 


A  -i  i  ^ 
ilE^f  (2)-  (CO)^ 


- T - X 

(n!)2  (a  +  CD)2 


CO 

r 


=■  I 
Inc  J  ^ 


(x)  exp 


ii) 


dX 


(VII- 65) 


where 


2  2CD 


a  +  CD 


(VII- 66) 


Denote  the  bracketed  term  by  G  (o): 

n 


G  (a)  =  - 

^  2nO 


j  H_^(x)  exp  (■  ^) 


dx 


(VII- 67) 


Substituting  into  Equation  VII- 45,  we  have: 


.  2E  \ 
r  ; 


o  N 


1  - 


o 


2a^  (CD)^ 
(a  +  CD) 


G“(a)C 
n  r 


n=0 


n! 


(\lI-08) 


Using  the  structure  of  the  Hermite  functions, ^  we  can  find  an  expres¬ 
sion  for  Gji(o). 


/  See  Cramer,  Reference  20. 


96 


^Ttliur  ZD.lLittU.Jnf. 


i 

i 

1- 
V 

‘Jt 

I  First,  u’c  evaluotr  G.j{'')  for  n  -  0,  2.  *,(i  and  deduce  the  expression 

f  for  iJ  rbitrary  n. 

I  Herrnite  Functions 

I  ^ 

^  H  (x)  =  1 

f  o 

I  (x)  =  x^  -  1 

^  H  fx)  =  x*^  -  6x^  +  d 

i  4 

H,(x)  =  x^  -  I5x^  =  45x^  -  15  G^(')  =  15(J^  -  (VII -69) 

f  6  6 

t 

Looking  at  the  sequence,  we  see  that 


G  (J) 
n _ 

G  (•)  =  1 
o 

G^  (-)  =  o'-’  -  I 
G^(o)  =  3(o"  -  1)^ 


and 


3 

o  -  1 


CD  -  a 
CD  a 


(VII -70) 


(VII-71) 


Ctoserve  that  for  c  -  1, 


God)  =  I 


(VII-7_>) 


G  (1)  =  0  n?«0 

n 


which  corresponds  to  the  case,  k  "  I . 
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Next,  evaluate  Cj^. 


N 

n  o 

=  -"ir 

r 


(VII- 73) 


From  Equations  VII-57,  VII-62,  and  VII-63  we  have; 


2  i  E  {|z|^}E^  .  D 
__ 


(VII-74) 


C  -  D 
C  D 


(VII-75) 


Substituting  Equations  VII-74  and  VII-75  into  Equations  VII-73  and 
VII- 68,  we  have: 


C  = 
n 


2jE{lzl^}E^  'j 


/  2|Eilzl^]Et  i 

(  D(C  -  D)” 


(VII-76) 


no  2E  /N 
r  o 


5  5  - 

2a^  (CD)  (2m)!  2 

(a-e'cD)  ^  ,2^' 

m'O  (m!2m) 


(VII-77) 


where  we  have  set  ni  =  -  since  the  odd  terms  in  the  series  were  zero. 
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_ \ _ 

iEUzI^lE  / 

N  1  +  BL 

o 


(Vll-gO) 


From  the  tabulation,  we  observe  that  there  is  very  little  difference  be¬ 
tween  and  d^^.  There  are  several  reasons  for  this  result: 

1 .  The  most  important  reason  is  that  we  have  considered  only  the  case 
where  the  target  has  zero  velocity.  As  shown  in  Figure  27,  the  target  is  exactly 
at  the  peak  of  the  scattering  function. 

If  the  target  had  a  non-zero  velocity  or  were  displaced  in  range  from 
the  peak  of  the  scattering  function,  the  difference  between  the  optimum  arxl  con¬ 
ventional  receiver  would  be  larger.  This  is  because  the  optimum  filter  uses  its 
knowledge  of  the  reverberation  scattering  function  to  partially  “tune  out”  the 
reverberation.  For  the  case  we  considered,  the  largest  amount  of  reverberation 
was  in  the  same  range-Doppler  location  as  the  target.  Thus,  the  optimum  filter 
could  not  tune  out  the  reverberation  peak  without  also  tuning  out  the  target.  For 
targets  away  from  die  peak,  the  optimum  filter  can  use  its  knowledge  more  ef¬ 
fectively. 


^TtltuT  B.littlr.inr. 

_ 


2.  A  second  reason  is  that  the  particular  scattering  function  we  have 
chosen  is  smooth  in  both  directions.  One  of  the  advantages  of  a  non- stationary 
iTKxlel  is  that  it  uses  the  non-uniform  distribution  of  scatterers  along  the  path 
of  sound  wave  to  improve  its  detection  capability.  Intuitively,  one  would  think 
that  the  more  non-uniform  the  distribution  is,  the  more  useful  knowledge  of  it 
would  be.  Thus,  our  choice  of  a  smooth  scattering  function  tends  to  negate  the 
effect  of  an  optimum  filter. 
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FIGURE  27  TARGET  LOCATION  WITH  RESPECT  TO  REVERBERATION 
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A 

a 

a(t) 


B 


H„(x) 

‘o<-) 

k 


L 


M 

X 

q 


X 


GLOSSARY 

attenuation  of  target  return 
length  piirameter  of  Gaussian  pulse 

intensity  functions;  average  number  of  scatterers  per  unit 
time 

rms  Doppler  shift  (reverberation) 
frequency  parameter  of  Gaussian  pulse 
S/N  ratio;  conventional  receiver 
S/N  ratio;  optimum  receiver 
received  energy 
transmitted  energy 
complex  envelope  ( f  (t)  =  u  (t)  e^^  ) 

n^  Hermite  function 

N/fodified  Bessel  function;  first  kind;  order  zero 
amplitude 

rms  length  of  scattering  hinction 
characteristic  function  of  scatterer  Doppler  shift 

characteristic  function  of  scatterer  velocity  is  a  random 
variable,  x  is  conditioning  variable 

height  of  ambient  noise  spectral  density  (double -sided) 
actual  additive  noise 
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"  % 


nA(t) 

N^(t) 

n^(t) 


D 


P  (x) 

*^0 

Q(a.b) 


R(t) 

r(t) 


Re 

S^<t) 

S„<t) 

S^(ui) 

Sj.(C) 

siv^.y^) 

t 

n 

u(t) 


Z 

n 

8 


Y 

A 


complex  envelope  of  siddltive  noise 

actual  reverberation  return 

complex  envelope  of  reverberation  return 

probability  of  detection 

probability  of  false  alarm 

probable  density  of  scatterer  Doppler  shift 

Marcum’s  Q  function 

actual  returned  signal 

complex  envelope  of  returned  signal 

real  part 

complex  envelope  of  desired  signal 

Fourier  transform  of  S_,(t) 

d 

complex  envelope  returned  from  n*^^  scatterer 
Fourier  transform  of  R(t) 
transmitted  signal 

characteristic  function  of  joint  scattering  density 
delay  due  to  n^^  scatterer 
actual  envelope 

complex  number  which  is  the  magnitude  and  phase  of  return 
phase  shift  of  target  return 

scatterers  per  unit  time  (uniform  case) 

a^  •»-  b^ 
a 
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i 

I 

T 


‘b 

X 

0(t) 

«„(t) 

e(T,u)) 


likelihood  ratio 


a  value  of  likelihood  ratio 
eigenvalue 

phase  of  transmitted  signal 

th  .  - 

n  eigenfunction 

two-dimensional  correlation  factor 
target  range  (delay  in  signal  return) 
carrier  frequency 
target  Doppler  shift 
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